Self-Consistent Renormalization Model of Mott Gap Collapse in the Cuprates 



R.S. Markiewicz 

Physics Department, Northeastern University, Boston MA 02115, USA 



A generalized antiferromagnetic approach to the Mott 
transition is analyzed with special emphasis on electron 
doped cuprates, where evidence for electronic phase separa- 
tion is weak or absent. Fluctuations are incorporated via 
a self-consistent renormalization, thereby deriving a 'nearly- 
antiferromagnetic Fermi liquid' susceptibility. The calcula- 
tion is sensitive to hot-spot effects. Near optimal doping, an 
approximately electron-hole symmetrical Mott gap collapse 
is found (quantum critical points). The calculation satisfies 
the Mermin- Wagner theorem (Neel transition at T = only - 
unless interlayer coupling effects are included), and the mean- 
field gap and transition temperature are replaced by pseudo- 
gap and onset temperature. The resulting susceptibility is 
used to calculate the doping dependence of the photoemis- 
sion dispersion, in excellent agreement with experiment. Dis- 
cussions of interlayer coupling, doping dependence of U, and 
extension to a three-band model are included. 



I. INTRODUCTION 

Schrieffer, Wen, and Zhang^ originally proposed that 
the magnetic insulating phase in underdoped cuprates 
could be understood via a spin density wave (SDW) 
approach to the Mott transition, and successfully de- 
scribed the spin wave spectrum of the undoped parent 
compound, which is an antiferromagnetic (AFM) insula- 
tor. Kampf and Schrieffer^ showed that precursors of 
the Mott transition could give rise to a pseudogap in 
the quasiparticle spectrum, between incipient upper and 
lower Hubbard bands (U/LHBs). Attempts were quickly 
made to go beyond mean field theories by incorporating 
fluctuation effects, but a number of problems soon arose. 
While some calculations found evidence for pseudogaps'^, 
others did not*. Many calculations found evidence for 
instabilities - either to incommensurate magnetism^ or 
to phase separation^. 

In a broader context, there is an ongoing controversy 
as to whether the Mott transition can be described start- 
ing from a band structure picture, or whether the strong 
correlations require a local point of view. For large hole 
dopings, a local picture seems justified by the observation 
that the magnetic correlation length ^ remains finite in 
the presence of the pseudogap, whereas mode coupling 
theories would predict a divergence in ^ as temperature 
T — > 0. However, this effect may be due to nanoscale 
phase separation physics. 

Newer experiments have suggested that indeed phase 
separation and/or stripe physics is present in the hole 
doped cuprates'''* down to arbitrarily small dopings^, 
thereby validating the early models. It remains difficult 



to develop a theory which simultaneously treats both the 
effects of strong fluctuations and (nanoscale) phase sepa- 
ration. A simpler alternative has recently been proposed. 
While phase separation is a significant complication For 
hole doping, this instability is greatly reduced or absent 
in electron- doped materials^°'^^, allowing a much simpler 
analysis. Moreover, for electron doping, the band picture 
involving short-range AFM order seems justified, in that 
the correlation length diverges for all dopings up to the 
QPT. 

Remarkably, simple mean field calculations suggest 
an explanation for both the electron-hole asymmetry^^ 
and the anomalous properties of the electron-doped 
cuprates^^'^"^. Thus, for a band structure with electron- 
hole symmetry, the AFM state is stable only at half fill- 
ing, being unstable against phase separation (negative 
compressibility) for any finite doping (this is also found 
in dynamic mean field theory^*). Introducing electron- 
hole asymmetry (here via a finite second-neighbor hop- 
ping t'), the doped system remains unstable for any fi- 
nite doping toward the Van Hove singularity (VHS) (hole- 
doping for t' < 0). However, for doping away from the 
VHS, electronic phase separation can be eliminated if 
the asymmetry is large enough. Similarly, the RPA spin 
wave spectra^ are found to be stable in electron doped 
materials, even when both upper and lower Hubbard 
bands cross the Fermi level^'^, while they are unstable 
against incommensurate SDWs in the presence of hole 
doping^^'^^. 

The resulting physics is much simpler for electron 
doping, with the gap decreasing gradually as the UHB 
fills. A quantum critical point (QCP) (Mott gap clos- 
ing) is found^° just beyond optimal doping, which can 
be understood theoretically^^ if the effective Hubbard 
U parameter decreases with doping (e.g., due to screen- 
ing). Since the pseudogap is associated with short-range 
AFM order, it appears predominantly at the hot spots - 
the points where the Fermi surface crosses the Brillouin 
zone boundary for long-range AFM order. As the gap 
shrinks, both upper and lower Hubbard bands cross the 
Fermi level, leading to two-band conduction, as observed 
exp er iment ally ^ ^ . 

The pseudogap closing also leads to a QCP near op- 
timal doping in hole doped cuprates, consistent with 
predictions^*'^^. However, its possible observation^^ is 
complicated by the presence of stripes, and this QCP 
has been analyzed in terms of a mode-coupling theory of 
two dimensional CDWs, taken as a model for stripes^^'^'^. 
Remarkably, for both electron and hole doping, the QCP 
occurs close to the point where hot spots vanish. It will 
be shown why this is expected to be the case^". 
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Thus, there appear to be alternative pathways to Mott 
gap collapse in the cuprates. Since the transition is bet- 
ter behaved in electron-doped cupratcs, it makes sense 
to develop a model which can describe the transition 
for these materials, then use the acquired insights to 
tackle the harder problem of hole doping with the atten- 
dant (nanoscale) phase separation. Here the mean field 
and RPA results are extended by incorporating fluctua- 
tions via mode-coupling theory^^, following Moriya's self- 
consistent renormalization (SCR)^^'^^ procedure. Mode 
coupling theories have been applied to charge density 
wave (CDW) systems^^'^^, and have led to a success- 
ful theory of weak itinerant magnetic systems^^'^^. They 
have also been used to study glass transitions^*, and re- 
cently extended to glasses in cuprates^^. The mode cou- 
pling analysis is particularly convenient, being the sim- 
plest model for which the Mermin- Wagner theorem is 
satisfied, and one can try to understand how to have a 
Mott (pseudo)gap without a superlattice. The present 
calculations are in general consistent with the results of 
Ref 3. The resulting pseudogaps compare well with re- 
cent photoemission experiments, while parallel results for 
hole doping provide clear evidence for additional compli- 
cations associated with stripe phases. 

The SCR theory involves a set of parameters^'' which 
also arise in the nearly antiferromagnetic Fermi liquid 
(NAFL)^'' and spin fermion^^'^^ models and in renor- 
malization group (RG) calculations of quantum phase 
transitions'^'' These parameters are generally esti- 
mated empirically from experiments. However, the good 
agreement between experiment and mean field theory for 
electron doped cuprates encourages us to try to calcu- 
late these parameters from first principles, in terms of a 
renormalized Hubbard parameter Uef / and a mode cou- 
pling parameter u. Most of the defining integrals are 
dominated by the region of hot spots. 

Some of these results have been reported previously in 
the discussion of the mean-field results^^'^'^ and in a con- 
ference procedings'^'"'. This paper is organized as follows. 
Section H presents an outline of the calculation, including 
the reduction to one band (Appendix A), doping depen- 
dence of U (Appendix B), and path integral formalism 
(Appendix C). The fiuctuation-induced correction to the 
susceptibilities is calculated, and is found to diverge for 
a two-dimensional system, driving the Neel temperature 
Tjv to zero (Mermin- Wagner theorem). Since the transi- 
tion occurs when a Stoner factor equals unity, it is con- 
trolled by the real part of the bare susceptibility. Hence 
Section HI reviews the properties of i?ex, showing that 
plateaus in % as a function of doping, q. or uj are all con- 
trolled by the physics of hot spots. In turn, these plateaus 
provide natural phase boundaries for QCPs^'^. The result- 
ing susceptibility has a form similar to that postulated for 
a nearly antiferromagnetic Fermi liquid (NAFL), and a 
calculation of the NAFL parameters (Appendix D) finds 
that there are extra (cutoff) parameters, which cannot be 
neglected. In Section IV, this renormalized susceptibil- 
ity is incorporated into the lowest-order correction to the 



electronic self energy, allowing a calculation of the spec- 
tral function associated with the pseudogap (Tjv = 0). 
Excellent agreement is found with the ARPES spectra 
of NCCO. Section V shows that inclusion of interlayer 
hopping leads to a finite Tjv (Appendix E). Results are 
discussed in Section VI, and Conclusions in Section VII. 

II. OUTLINE OF THE CALCULATION 
A. t — J vs Mode-Coupling t—U Models 

While the strongly correlated Hubbard model is often 
approximated by a t — J model, this is not appropriate 
in the present analysis. In the electron doped cuprates, 
both Hubbard bands need to be accounted for, since (a) 
ARPES can detect both bands (at least up to the Fermi 
level), and (b) the Mott gap is found to collapse with 
doping, leading to an overlapping of both bands with the 
Fermi level. This situation is difficult to incorporate into 
the t — J model, where one Hubbard band is neglected. 
For this reason, an alternative procedure is used in the 
present paper. 

While mean field calculations reproduce the low tem- 
perature properties, they do not account for thermal fluc- 
tuations, and lienc;e predict that Neel order persists to 
too high temperatures - indeed, for a two-dimensional 
system, Neel order can exist only at T = (Mermin- 
Wagner theorem)''®. Even when fluctuation effects de- 
press the onset of long-range order, short-range order 
can still persist up to the mean-field transition tempera- 
ture, producing a pseudogap similar to the ones found in 
one-dimensional CDWs^^. In the AFM these fluctuations 
are associated with spin wave excitations, in particular 
the Goldstone modes near Q. Here, this effect is calcu- 
lated, following the self consistent renormalization (SCR) 
scheme of Moriya^^'^^. 

The calculation recovers the mean- field results, that 
antiferromagnetism is stable for electron doping (in the 
presence of a second- neighbor hopping parameter t' < 0), 
but for hole doping the lowest energy state is either an 
incommensurate SDW^^ or phase separated^^'^^. [Unre- 
stricted Hartree-Fock calculations find that for t' = Q, 
the incommensurate antiferromagnetic state has lower 
free energy, while for sufficiently negative t', the phase 
separated state is the lower free energy state^°.] While 
the SCR technique can be generalized to deal with com- 
peting phases'*^, only the antiferromagnetic fluctuations 
will be treated here. To keep the results consistent with 
a low-doping t—J analysis, the mode coupling parameter 
is chosen to reproduce the t — J expression for the spin 
stiffness parameter at half filling. 
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B. Model Dispersion and Doping Dependence of U 

The cuprates are treated in a one-band model. By 

comparison with a 3-band model (Appendix A) , this can 
be shown to be an excellent approximation for the mag- 
netic properties. The bare electronic dispersion is 
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with Cj = cos kiU. The dispersions for undoped 
Sr2Cu02Cl2 SCOC and electron-doped NCCO can be fit 
by assuming t = 0.326eV, t' /t = —0.276, with U taken as 
an effective doping dependent parameter^^, with U — 6t 
at half filling. Remarkable, virtually the same parame- 
ters are found''^ to describe the spin wave spectrum'*^ in 
La2Cu04: t = OMeV, t' /t = -0.25, and U/t = 6.2. The 
former values will be used here. 

Many textbooks on strong correlation physics^**'^^ note 
that the Hubbard U should be doping dependent, based 
on the original results of Kanamori'*^, but there arc no 
satisfactory results for the doping dependence in the 
cuprates. A simple model calculation, which gives semi- 
quantitative agreement with experiment in NCCO^'^'^"'^, 
is described in Appendix B. 



C. Mode Coupling Calculation: 
Derivation of NAFL Parameters 

A naive perturbation theory breaks down due to mode 
coupling, and the SCR scheme is introduced to calcu- 
late the renormalized susceptibility near the antiferro- 
magnetic wave vector Q. The (path integral) formal- 
ism is reviewed in Appendix C, and only the main re- 
sults are given here. The divergence of the suscepti- 
bility is controlled by the (inverse) Stoner factor Sq = 
1 — Uxo{Q + Qt^), where xo is the bare magnetic sus- 
ceptibility. Within the SCR, So is renormalized to 5 > 0, 
so there is no divergence, but when ^ = 0+ there are 
strong fluctuations (pscudogap regime). The physics is 
controlled by the dispersion of 6 near Q, 



5q{ij) = 5 + Aq^ + A^ql - Bo? - iCuj, 
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Eq. C18; the important parameters A, B, and C are 
evaluated in Appendix D (B is small and can generally 
be ignored). This leads to a susceptibility 



X(9,'^) = 



XQ 



as a possible mechanism for generating a finite Ncel tem- 
perature. The similarity of Eq. 3 to the corresponding 
result for CDW's^^ should be noted - the SCR is a form 
of mode coupling theory. 

Equation 3 is of nearly antiferromagnetic Fermi liq- 
uid (NAFL)^" form, and the same parameters also arise 
in a renormalization group calculation of QCP's"^^'^^. 
Thus, evaluation of the parameters A, B, and C would 
amount to a microscopic derivation of NAFL theory. 
The present work evaluates these parameters in terms of 
two interaction parameters, the Hubbard U and a mode 
coupling parameter u. It is found that U has an im- 
portant doping dependence, estimated in Appendix B, 
which is consistent with experiment. However, an at- 
tempt to directly calculate the mode coupling parame- 
ter u fails, giving anomalously small values - this prob- 
lem has been noted previously, although there is debate 
about whether u diverges^^ or vanishes^^'^''. Here, the 
t — J model is used to fix the value of u. Finally, it 
must be noted that the SCR theory is a model of weak 
ferro- or antiferromagnetism^^ , and is here found to un- 
derestimate the Mott gap splitting near half filling. This 
most likely is due to imperfect self-consistency: the pa- 
rameters arc estimated from the bare susceptibility 
while the opening of the pseudogap leads to significant 
modifications of %• 

Despite these limitations, the resulting model is well 
behaved, and in good qualitative agreement with exper- 
iment, suggesting that a full derivation can ultimately 
be carried out along these lines. There are a number of 
deviations from the simplest NAFL theory. In Section 
HI, it will be shown that the theory contains two ad- 
ditional CMtoj(f parameters, Qc and w", which cannot be 
neglected (or sent to infinity). One consequence of this is 
that the A parameter develops a significant temperature 
dependence, particularly in the electron-doping regime. 

Equation 3, with Imx ~ i-^, only holds in the presence 
of hot spots'^^. Hot spots complicate the evaluation of 
the A, B, and C parameters, since the expressions contain 
integrals which are formally divergent near the hot spots. 
For the present band structure, hot spots exist only when 
the chemical potential /i is in the range 4t' < /i < 0, 
or for doping 0.25 > a; > —0.19 (electron dopings are 
considered as negative). 

The following section demonstrates how hot spots con- 
trol the properties of the magnetic susceptibility (plateau 
formation), and summarizes the evaluation of the SCR 
parameters. The integrals arc more singular at the end 
jpoints of the hot spot range - the H and C points. 
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with coefficients given by Eqs. C39-C42 in terms of A, 
B, and C. Here and below q and Q arc treated as two- 
dimensional vectors, while q^ is introduced explicitly, and 
ttz = Az/A. Interlayer coupling will be ignored, az = 0, 
until Section V and Appendix E, where it is considered 
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III. HOT SPOT PLATEAUS AND GENERIC 
QCPS 

A. Plateaus in Doping Dependence 

1. The Pseudo-VHS 

The present analysis is based on a self-consistent renor- 
malization scheme: in two dimensions, fluctuations pre- 
vent the establishment of long range AFM order. Hence, 
the relevant quantity on which the study is based is the 
bare magnetic susceptibility. 



Xo 



(4) 



where (5 is a positive infinitesimal. This susceptibility has 
been analyzed in a number of papers, but generally only 
Im{x) is explored in detail (e.g., Refs. 47-49), whereas 
the Stoner criterion involves Re{x), which was studied 
in Ref. 50. The extended discussion which follows is 
intended to bring out salient features for the computa- 
tion of the NAFL parameters. The doping dependence 
of XoiQ:'-^) is illustrated in Fig. la, where Q = {n,Tr). 
The susceptibility has a remarkable doping dependence, 
with the large peak at the Van Hove singularity (VHS) 
shifting^'' to half filling with increasing temperature T. 
The peak position of this 'pseudo-VHS' defines a temper- 
ature Tv{x), Fig. Id (circles). This behavior can readily 
be understood from the form of XoiQ^O)! Eq. 4. The de- 
nominator = —^t{cx + Cy), is independent oft', 
and hence has a stronger divergence than the density of 
states (dos). Indeed, this divergence matches the strong 
VHS found for t' — (perfect nesting), and like that VHS 
falls at half filling, a; = 0. There is one crucial difference 
- at low temperatures, this divergence is cut off by the 
Fermi functions, which leave the integrand non zero in 
a wedge which intercepts the zone diagonal (where the 
denominator vanishes) only at isolated points: the hot 
spots. Hence, the residual divergence at low T is still 
dominated by the conventional VHS. However, at finite 
T, excitations along the zone diagonal become allowed, 
leading to a stronger divergence of xo('3, 0) near x — 0. 

The strong temperature dependence of the pseudo- 
VHS is in strong contrast to the density of states, Np, 
Fig. lb, and also with the pairing correlations^". The de- 
nominator of the pairing susceptibility involves the sum 
of the energies, + £j:^Q = —&t'cxCy, rather than their 
difference (as in Eq. 4), and hence always peaks at the 
ordinary VHS. 

The difference between nesting and pairing susceptibil- 
ities has a fundamental significance. By mixing electron 
and hole-like excitations, the superconducting gap is al- 
ways pinned to the Fermi level, and can open up a full 
gap at any doping. On the other hand, a nesting gap 
need not be centered on the Fermi surface, and is con- 
strained to obey Luttinger's theorem, conserving the net 



number of carriers in the resultant Fermi surface. Hence, 
the only way a nesting instability (such as antiferromag- 
netism) can open a full gap at the Fermi level is for the 
instability to migrate with increased coupling strength 
to integer filling of a superlattice zone (e.g., half filling of 
the normal state). 




FIG. 1. (a) Susceptibility xo at Q as a function of doping 
for several temperatures. From highest to lowest curves near 
X = 0.1, the temperatures are T = 1, 100, 300, 600, 1000, 
2000, and 4000 K. Dotted line = 1/Ueff, dot-dashed line = 
1.5/f/e// • (b) Density of states Nf for the same temperatures, 
(c) Susceptibility xo at Q as a function of doping for several 
frequencies a.t T = IK: u = 0.01, 0.1, 0.3, 0.6, 1.0 eV. (d) 
Pseudo-VHS (peak of xo) as a function of temperature Tv 
(circles) or scaled frequency = cj^/tt (squares); triangles 
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Since the susceptibility has such a distinct temperature 
dependence from the density of states, one might ask how 
the frequency dependence compares. This is illustrated 
in Fig. Ic at low temperature (IK). While the frequency 
introduces additional sharp features and has an overall 
very distinct appearence from the T-dependence, never- 
theless the main peak also shifts from the VHS toward 
lower doping with increasing a; - in fact, the shift is al- 
most the same when comparing fiuj and nkBT, Fig. Id. 
The dashed line in Fig. Id is T~ = huj~ /irkB, with^^ 



4t(A - r) 



1 



(5) 



with r = 2t'/t and jj, — /i/2i. The proportionality of 
frequency and temperature dependences holds only in the 
hole doped regime: temperature shifts the susceptibility 
peak only to half filling, a; = 0, while frequency will shift 
the peak beyond half filling {x < 0). 

Also in contrast to Np, the susceptibility has (at low 
T) a plateau shape, with sharp falloff in intensity beyond 
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the plateau edges on both electron and hole doping sides 
of half filling. This shape is characteristic of hot spot 
physics. Hot spots are those points where the Fermi sur- 
face (FS) intersects the replica FS shifted by Q. They 
are located at c^; = — = Cxo, with 

Cxo = cos afc^o = \J^' 

and equivalent points. The edges of the plateau are those 
points at which the overlap terminates (hot spots cease 
to exist). For the band structure of Eq. 1 these points 
occur at chemical potential fj, = At' = /ii, (the VHS) 
and 0, or at dopings x = 0.25, -0.19 (taking electron 
dopings as negative). Since these two end points play 
an important role, it is convenient to label them, and 
they are here called 'hot' hot spot and 'cold' hot spot 
(or H-point and C-point) for the hole and electron-doped 
termination points, respectively. It will be demonstrated 
below that at each doping, the hot spots also lead to 
a susceptibility plateau in momentum space, around Q, 
collapsing to a logarithmic (square root) divergence at 
the H- (C-)point. The iJ-point is the VHS, and hence 
also involves a conventional ETT. The physics is simpler 
near the C-point, where the topology hardly changes but 
the FS and Q-FS become decoupled (it is therefore a form 
of Kohn anomaly^°). 

2. Mean Field Mott Transition 

For the parameter values expected in the cuprates, 
these susceptibility plateaus control the physics of the 
Mott gap collapse^°. As a function of doping, the mean 
field Mott gap is found to close at a doping just beyond 
the edge of the plateau, for both electron and hole dop- 
ing, Fig. 2. The solid and long dashed lines are the com- 
mensurate and incommensurate mean field Mott transi- 
tion temperatures T*{x) calculated using the estimated 
Ueff{x), dotted line in Fig. 1. For electron doping, there 
is a double transition, first from commensurate to incom- 
mensurate antiferromagnetic order at the plateau edge, 
then to the loss of any magnetic order at a slightly higher 
doping (inset a). For hole doping, the dominant antifer- 
romagnetic order is incommensurate for all dopings, but 
the difference in T/v becomes significant only near the H- 
point (inset b). When fluctuations are included (below), 
it is found that the Neel transition is shifted to zero tem- 
perature, while a pseudogap first appears near the mean 
field Tjv. Interlayer coupling can then restore a finite Tjv, 
Section V. For the real cuprates, the terminations of the 
Mott gaps are preempted by superconducting transitions, 
close to the critical regime. 

While the present local calculation (minimizing the 
free energy at a fixed doping) finds incommensurate mag- 
netic order for hole doping, some global calculations 
(comparing free energies over a wider doping range) find 
that this is precluded by a phase separation instability. 
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FIG. 2. Mean field magnetic transition temperatures de- 
termined from Stoner criterion using Ueff of Fig. 1. Solid 
line: commensurate (at Q); long dashed line: incommensu- 
rate. Dot-dashed fine = lOTjv, where Tjv is the onset of long 
range AFM order, from [ 51] and [ 52] (with filled circles). 
Insets = blowups near C- and H-points. Squares in inset b — 
pseudogap data of [ 53]. 

The structure in the low temperature susceptibility. 
Fig. 1, with its largest peak at the H-point on the hole 
doped side, is in striking contrast to the calculated dop- 
ing dependence of the Neel transition. Fig. 2, which has 
a broad plateau on the electron-doped side, but falls 
off more quickly with hole doping, showing no sign of 
a peak near the VHS. This contrast can be accounted 
for by two effects. First, the shift of spectral weight 
with temperature of the pseudo-VHS, noted in Fig. 1, 
would tend to produce a symmetric falloff of Tpf with 
either electron or hole doping. But the dos peak at 
the VHS leads to better screening of Ueff for hole dop- 
ing, thereby further depressing T^. The experimentally 
observed^^ T/v (dot-dashed line) shows an even stronger 
falloff with hole doping, perhaps due to phase separa- 
tion. Since stripes can frustrate magnetic order, the fig- 
ure also includes the magnetic ordering temperature of 
quasi-static stripe arrays, from Nd-substituted LSCO^^, 
which is taken as a lower bound for the Neel ordering 
transition in the absence of stripes. The mean field cal- 
culation provides an approximate envelope of the result- 
ing data, but overestimates the transition temperatures 
by a factor of 10. Note that in the hole doped regime, 
there is good agreement between the mean field transi- 
tion and the pseudogap^° (squares in Fig. lb — data of 
Krasnov^'^, assuming 2A = 4.6T*). Calculation of the 
Neel transition beyond the mean field level will be dis- 
cussed in Section V. 
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B. Plateaus in Momentum Space 
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1 . Plateaus 
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FIG. 3. Susceptibility xo near Q for a variety of dopings at 
T = lOOA'. From highest to lowest solid curves near S = Q, 
the chemical potentials are jj, = -0.35, -0.30, -0.25, -0.20, -0.15, 
-0.10, -0.055, -0.02, and eV. For the dashed curves (top to 
bottom), /i = -0.352, -0.355, and -0.359eV. 

In analyzing either thermal fluctuations or the quan- 
tum fluctuations associated with QCPs, it is necessary to 
understand the susceptibility near the AFM vector Q. At 
each doping, hot spot physics leads to a plateau in mo- 
mentum space, centered on Q. Figure 3 shows how xo 
varies near Q at a low temperature (lOOK) for a series 
of different dopings. Results near T = are presented in 
Ref. 35. For all dopings there is a plateau in q. The width 
of the plateau at T = can be readily determined^": in 
any direction, it is the minimum q needed to shift the 
replica FS so that the hot spots are eliminated. This 
can be found from the dispersion, Eq. 1, by substituting 
k->{Q + q)/2, or 

-2t{Sx + Sy) - At'SxSy = fl, (7) 

with Si — sin(gia/2). As shown in Fig. 4, this for- 
mula agrees with the (anisotropic) plateau width mea- 
sured from Fig. 3 (circles). The inset shows the shape 
of the plateau as a function of doping. The diamond 
shape of the plateau, Eq. 7, is related to the profile of 
the hole pockets formed by the overlap of the shifted and 
unshifted FSs. Specifically, the plateau is the region of 
overlap of the two hole pockets, shifted to have a common 
center, as illustrated in Fig. 5. The remaining parts of 
the pockets also show up, as ridges^'' in the susceptibility, 
radiating from the corners of the diamond (similar to the 
peaks in the /i — O.ObeV data in Fig. 6, below). As noted 
by Benard, et al.^^, the susceptibility in two-dimensions 
acts as a FS caliper. The plateau width leads to a natu- 
ral limit on the magnetic correlation length, l/^c in 
agreement with experimental data^^ (squares in Fig. 4), 
as noted previously^^^'^". 
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FIG. 4. Plateau width Qc, comparing Eq.7 (solid lines) and 
the measured widths (circles) from Fig. 3. Upper curve along 
[qc,0] direction, lower along [qc,qc]/V2 direction. Squares = 
experimental inverse correlation lengths from Ref. 55; dia- 
monds = TX/5000K . Inset — plateau boundary for a series of 
chemical potentials n from (smallest) to -0.359eV (largest). 
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FIG. 5. Illustrating origin of plateaus (dotted line) from 
crossed hole pockets (short dashed lines). 

The plateaus in q help in understanding the doping de- 
pendence of the susceptibility near Q, Fig. la. At each 
doping on the plateau (in x) there is a plateau in q cen- 
tered at Q, with the width of the plateau decreasing to 
zero as a; ^ xc, Fig. 6. The critical points q = qc are 
precisely those points at which the Q + (f-shifted-FS no 
longer overlaps the original FS. 
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FIG. 6. Susceptibility xo near Q for several dopings near 
the C-point. Upper group a.t ^ = -0.05eV, middle at ^ 
= (C-point), and bottom aX ^ = +0.05eV. Temperatures 
are T = (dotted lines), lOOK (short dashed lines), 

lOK (long dashed lines), IK (solid lines). Horizontal line = 
(7.//(m = 0). 
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FIG. 7. Bare susceptibility xo, for a; = 0, (beaded line); 
various solid lines show contributions of individual quadrants, 
shifted vertically by 0.1 for clarity. The special points in 
the Brillouin zone are F = (0,0), X = (7r,0), and S — Q. 
(T=100K.) 



2. Cusps 

For electron-doping beyond the C-point (/i > 0), the 
plateau ends and the susceptibility displays split peaks 
away from Q, Fig. 6, with a dip in between. The change 
in character of xo nieans that ii ~ is a QCP. (The cor- 
responding QCP at the H-point was analyzed in Ref. 50.) 
However, the magnitude of xo also changes rapidly near 
/i = 0, so there should be a transition to a non-magnetic 
phase near the same doping^", as discussed in the pre- 
vious subsection (note the Hne depicting U{fi = 0) in 
Fig. 6). 

The origin of these /i > cusps can be readily under- 
stood from Fig. 7. Here, the contributions of individual 
quadrants to the x integral are plotted separately, with 
each quadrant containing two hot spots. It can be seen 
(Appendix D) that the dispersion at each hot spot con- 
tains a cusp, near which the dispersion is linear in \q\. 
However, when adding the contributions of the 8 hot 
spots, the linear terms cancel, leaving a quadratic dis- 
persion. For n > and q — there are no hot spots - 
the Fermi surfaces in k and k + Q do not intersect. How- 
ever, translating one Fermi surface by q will lead to an 
intersection, with corresponding hot spot, beyond some 
threshold qc- Since only one or two hot spots are restored 
for a given direction of qc, the linear terms do not cancel, 
leading to a linear in q dispersion for q > qc- Exactly at 
M = 0, Xo has a ^/q cutoff as q Q (Appendix D3). 



Technically, similar cusps also arise at the plateau 
edges^° for electron doping, > ^ > —0.22eV. The tops 
of the plateaus are not completely flat, Fig. 8a: the high- 
est susceptibility is shifted away from Q, and the sharp 
steps near Q are again hot spot effects, this time associ- 
ated with the loss of hot spots at large values of q — Q. 
However, these effects are much weaker than those asso- 
ciated with /X > (Ax/x ^ 0.5% - compare the vertical 
scales of Figs. 6, 8). Thus near the mean-field transition 
any structure on the plateaus is smeared out by thermal 
broadening. Even at T = 0, these features are likely to 
be negligible compared to dispersion in U which arises 
from renormalization effects^^. 



1.228 




FIG. 8. (a): Expanded view of susceptibility xo on the 
plateaus near Q for a variety of dopings at T = lOOA" (solid 
curves^ or IK (dashed curves). From highest to lowest curves 
near Q, the chemical potentials are fi — -0.20, -0.15, and -0.05 
eV (for both solid and dashed curves). All curves except 
fi = — 0.20eV^ have been shifted vertically to fit within the 
expanded frame, (b): Similar plateaus for the hole doped 
materials (T = IK), with (from highest to lowest) /i — -0.359, 
-0.35, -0.3, -0.25, and -0.22 eV. 
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3. Curvature (A) 

The plateau is a region of anomalously small local cur- 
vature A = A/U (Eq. 2) of the susceptibility, XoiQ + q) = 
XQ — Aq^, where A is an important NAFL parameter. 
Clearly, at T = lOOK the curvature A has gone neg- 
ative near the iJ-point, Fig. 3. At even lower tem- 
peratures, it reverts to positive values, Fig. 8b. The 
temperature dependence of the normalized parameter 
A' = {TT/aY{U /t)A is illustrated in Fig. 9 at several dop- 
ings. The temperature dependence is dominated by di- 
vergences at both H- and C-points. The divergence at 
the H-point Fig. 9a is the well-known logarithmic VHS. 
However, at finite temperatures spectral weight is shifted 
away from the VHS and A turns negative, only recovering 
a positive sign above T ~ 2000-?^. The temperature at 
which A turns negative can be defined as Ti„comm- ^ < 
for T > T^ncomm- From Fig. Id, Ti„comm is comparable 
to but larger than Ty (for x < 0.06 A remains posi- 
tive). This in fact explains the origin of Tincomm- Fig- 
ure 9a demonstrates that A is negative at T ^ beyond 
the H-point {/Li = —OAeV). Thus, increasing T above 
Ty produces the same susceptibility crossover. A similar 
crossover was discussed by Sachdev, et al.^^, except that 
they assumed that in the high temperature phase the 
AFM fluctuations remained centered on the commensu- 
rate Q, whereas here A is negative. At sufficiently high 
temperatures A again becomes positive for all dopings - 
i.e., the leading singularity of xo is always at Q. 

At the C-point, the collapse of the plateau width trans- 
lates into a divergence of the curvature at Q (A — > oo). 
This divergence of the high-temperature susceptibility is 
cut off at low T, Fig. 9d, when the thermal smearing be- 
comes smaller than the plateau width. For smller T, A 
is controlled by the curvature on the plateau. The tem- 
perature at which A has a peak, defined as T^, is plotted 
as diamonds in Fig, 4 (the peak is only found for a; < 0). 
Rather surprisingly, scales with the plateau width qc, 
even though the dynamic exponent is z = 2. Further, the 



1.5 



maximum slope scales approximately as A„iax '~ -^a 
which follows from the fact that A ~ q_ 
point. 

At intermediate doping. Fig. 9b, c, A is generally a 
scaled-down version of the behavior near the two end 
points, with a crossover near fi = — 0.25ey, where the 
T-dependence is weak. Also for intermediate tempera- 
tures, there can be fine structure on the plateau (e.g., 
solid lines in Fig. 8a) which can lead to wild swings in 
A{T). However, at these dopings they are not relevant, 
since the susceptibility peaks are away from Q, and this 
fine structure is not generally reported in Fig. 9. 




12 3 

log,„(T(K)) 



12 3 4 

log„(T (K)) 



FIG. 9. Temperature dependence of A' for several dopings. 



C. Plateaus in Frequency 
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FIG. 10. (a) /mx(Q,a;), (b) Rex{Q,t^), (c) Jmx/w = C, 
and (d) dRex(Q , ^) / duj , for (a,c): /i = (solid line), -0.05 
(long dashed line), -0.10 (dashed line), -0.15 (dotted line), 
-0.20 (dot-dashed line), -0.25 (dot-dot-dashed line), and 
-O.SOeV (short dashed line); (b,d): a; = (sohd line), 0.04 
(long dashed line), 0.10 (dashed line), and 0.15 (dotted line). 

Figure 10 illustrates Imx{Q,io), Rex{Q,oj), and 
Imx/i^ = C . At T = 0, the imaginary part of the sus- 
ceptibility x{Qi^) can be calculated analytically''^: 

/m(x(Q, io)) = Y.{f{e^) fiH+Qmn+Q H - ^) 

k 

_ F(0i,fc)-F(02,fc) 



At 



(8) 



where F{6,x) is an elliptic integral, k = -y/l — (uj/St)'^, 
and sin(Oi) = sin{(pi)/k, with 



cos (01 



if UJ < LO^ 

uj/2 if w > 



(9) 
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if LO < Uq 
if LU > LOq ' 



(10) 



and 



with jl = /i/2t, cD = uj/4:t, Cj_ = a± ^ y u.± 

a± = 1 — {p, ± Co)/t. The real part Rex can be found 
from the Kramers-Kronig result, 



IT 



Imx{Q, uj')u}'duj' 



(11) 



o 




FIG. 11. C calculated for several values of /x: jx = -0.355 
(diamonds), -0.357 (circles), -0.358 (squares), -0.359eV (trian- 
gles) [fiv = -0.3599eV]. Inset: Band dispersion eg (solid line) 
^k+Q (dashed line), for fi — 0. Arrow = ujc ■ 

Thus, there are also plateaus in the frequency depen- 
dence of Rex- Furthermore, hot spots generate an imag- 
inary part of the susceptibility linear in frequency, which 
also approximates a plateau, particularly near the H- 
point. Fig. 11. The origin of this plateau and of the 
critical frequencies uj~ , Eq. 5, and 



8t 



UJQ 



(12) 



can be understood from Fig. 12. The thick (thin) solid 
lines represent the original (Q-shifted) Fermi surfaces, 
while the dashed lines represent 



(13) 



for various values of uj. Equation 13 gives the points at 
which the denominator of xo(Qi Eq. 4, vanishes. Thus 
at T = 0, /m(xo(Qj'^)) is proportional to the length of 
the dashed line lying between the original and Q-shifted 
FSs (i.e., where /(eg)-/(eg_^-) = ±1). Since the two FSs 

meet at an angle, forming a wedge, Im{xo{QT^)) ~ 
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FIG. 12. Origins of critical cutoffs. Thick solid line = FS; 
thin solid line = Q-shifted FS; dashed lines — Eq. 13, for 
several values of ui. Chemical potential /i = (a) 0, (b) -0.1, 
(c) -0.14, (d) -0.2eV. horizontal arrows indicate u>o, vertical 
arrows lu~ . 

From Fig. 12, the critical frequencies (denoted by ar- 
rows) are points where the uj dependence of this length 
changes abruptly, leading to a sharp change in Imx- 
Thus, near the H-point, the plateau width is lo^ (in- 
set. Fig. 11), while near the C-point it is ujq. The ver- 
tical arrows in Fig. 12 indicate w", where the dashed 
line (Eq. 13) intersects the FS at the zone boundary, 
while the horizontal arrows are"*^ ujq, where the dashed 
line ceases to intersect the Q-shifted FS. There is a 
crossover at /ic — —O.lAeV: for /i > fic, wo < 
while for /i < fic, > ■ Combining Eqs. 5,12, 
ujo = oj- at ^c. = [l-z(2- Vz)^]2</t = -0.1384ey, with 
z — 1 — T. For UJ > min{uio,uJc}j -^"z(xo(Qj w)) ~ lj-^^^, 
so C ~ l/w^/^ - i.e., the susceptibility is no longer on 
the plateau. 

The height of the plateau C is an important parameter 
of the SCR model. It can be represented as another fre- 
quency wi = 1/C, with C = UC{uj = 0). From Eq. C22 
of Appendix C2, C can be found explicitly 



1 



2t:JsIq{1 + Tc^o) 



1 



(14) 



(with J = 4t^/U, si 



1 ~ c^o). Defining a width 



parameter a^j — mm{a^,a"}, with — ojq/uji, then 

_ 27rt(l-T) 1+TCrf) 1 

= 77 [ ] = — • (15) 

This latter is in good agreement with the numerical re- 
sults (arrows in Fig. 12) and is similar to the result found 
by Onufrieva and Pfeuty^*^, using a hyperbolic band ap- 
proximation valid near a VHS, loi/uj'^ — 27rt(l — t)/U . 

Because of the dynamic scaling lu q^, this crossover 
is also reflected in the behavior on the plateau in q, Fig. 8: 
for fi > — 0.14ey, the plateau has a negative curvature, 
which can almost be scaled between different dopings, 
while for fj, < — 0.14ey, the plateau starts to fill in, ul- 
timately developing a peak at Q. (See also Fig. 3a in 
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Ref. 35.) Note that the plateau width collapses in fre- 
quency at both the H- and C-points, while the collapse 
in wave number {qc 0) is only present near the C-point. 



FIG. 13. Calculated values of A (circles) and C (squares), 
for U = 6t. Solid line = Eq. 14 



D. Parameter Evaluation for Mode Coupling Theory 

The evaluation of the SCR parameters A and C was 
discussed above. The collapse of the q and/or uj plateau 
widths near the H- and C-points leads to the introduction 
of additional parameters qc and a^^ . The narrow width of 
the g-plateau, particularly for electron doping, leads to 
an additional complication not included in the conven- 
tional SCR analysis: the curvature of the bare suscepti- 
bility near Q = (tt, tt) (the S-point of the BZ) is strongly 
temperature dependent, and for some dopings may even 
change sign. In principle, it is not difficult to incorpo- 
rate an A{T) into the analysis near the mean-field Nccl 
temperature (pseudogap onset). But for the present 
2D system, long range Neel order only sets in at Tjv — 0, 
and for T << T^, a self consistent value of A should be 
found, by taking into account the effect of the pseudo- 
gap in modifying the electronic dispersion and hence x- 
For the present, this complication is ignored, and in the 
following section A is taken a,s A — A{T^), where 
is the magnetic pseudogap onset, the temperature where 
Xo{Q)Ueff = 1, using the effective Ueff found earlier-'^^ 
(Appendix B). This should be the most important A for 
controlling the pseudogap, and moreover at lower tem- 
peratures the band renormalization should strongly mod- 
ify A{T). With this choice, the resulting A{^) is plotted 
in Fig. 13, along with the C parameter, evaluated at 
T = 0. Note that for electron doping, this choice of A 
is always positive and varies smoothly with doping, di- 
verging at the C-point. By contrast, for hole doping A 
is often negative, again illustrating the instability of the 
uniform AFM phase. [In principle, a positive A can be 
found by taking an incommensurate nesting vector; for 
this paper, only commensurate nesting is considered; the 
negative A will suggest when electronic inhomogeneity 
may be important.] Given A and C, Fig. 14 shows the 
calculated values of Xq and ujsf, normalized to 

10 ■ , ■ , ■ , ■ , — ^ 
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FIG. 14. Calculated values of Xq/C^ (solid line) and ujgf^^ 
(short dashed line), assuming U — 6t. Long dashed line = 
doping a:: (/i) (xlO). 

In the following, the present results are applied to 
understanding the ARPES spectra of electron doped 
cuprates, concentrating on the four dopings analyzed by 
Armitage, et al.^". For convenience. Table I summarizes 
the parameters for these dopings. From the mean-field 
analyses^^, the effective Hubbard parameters were found 
to be Ueff/t = 6 (a; = 0), 5 (a; = -0.04), 3 (x = -0.10), 
and 2.5 {x — —0.15). [These numbers differ somewhat 
from those of Ref. 11, since a second neighbor hopping, 
t", was included in the latter analysis, to give the best 
fit of the Fermi surfaces.] The Stoner factor has a quan- 
tum correction, 77, Eq. C34, which tends to suppress the 
AFM transition; hence a smaller renormalization of U is 
required. This is reflected in Table 1: for x = -0.1, -0.15, 
there are two rows, the upper row using the mean-field 
U parameters, the lower with the quantum correction. 
Note that the [/'s are enhanced by essentially the quan- 
tum correction factor. These values will be used in the 
subsequent analysis. 

The SCR analysis also involves a mode coupling pa- 
rameter u, evaluated in Appendix D6. As found pre- 
viously, direct evaluation of this parameter is unsatisfac- 
tory - the results of Table 1 being anomalously small due 
to the flatness of the susceptibility plateau (dx/doj ~ 0), 
Fig. lOd. Below, an empirical way to estimate u is sug- 
gested. 



Table I: Electron Doped Cuprates 



X 


U/t 


A/a^ 


wi(eV) 




qca 


V 




(cV) 





6 


0.696 


0.345 


0.583 


0.635 


1.29 


1020 


760 


-0.04 


5 


1.16 


0.540 


0.455 


0.518 


1.25 


850 


3200 


-0.10 


3 


1.34 


1.32 


0.176 


0.342 


1.19 


500 


2700 




3.5 


1.56 


1.13 


0.206 




1.24 




2300 


-0.15 


2.5 


1.75 


2.16 


0.054 


0.172 


1.08 


56 


4000 


75 


2.9 


2.03 


1.86 


0.062 




1.15 


55 


3500 



It is convenient to compare the present results with 
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parameters estimated for the SCR model^^ from experi- 
mental data for (optimally) hole-doped cuprates. The pa- 
rameters are defined as Tq — Aqg/2'KC, Ta — Aq^ /2xo, 



yo = So{T = 0)/Aql, and yi ~ Ua^u/n^AC. The re- 
sults are listed in Table II, where the first line gives the 
hole-doped results estimated in Ref. 29. Moriya, et al.^^ 
took q% = l/Ana^ {qsa = 0.282), while for Table 11 it is 
assumed that qs = Qc- A key difference is that Moriya, 
et al.^^ assume the system is in the paramagnetic phase 
{yo > 0) at and above optimal (hole) doping, while in 
the present work y^ < 0, and the system is paramagnetic 
due to the Mermin- Wagner theorem, with the Mott gap 
appearing as a pseudogap. The small magnitude of yo is 
suggestive of a system pinned close to a QCP. Finally, the 
parameter yi is estimated using the value — 0.256eV 
(below), and not the anomalous values of Table 1. 

Table II: SCR Parameters 



X 


To (K) 


Ta (K) 


yo 


yi 


- 0.2 


1600-4000 


3000-10000 


0.01-0.02 


3 


0.0 


180 


1150 


-5.27 


0.75 


-0.04 


310 


1300 


-3.31 


0.7 


-0.10 


380 


670 


-1.23 


1.5 


-0.15 


200 


220 


-0.31 


1.85 



is presumably a deficiency of the present model in not 
using fully self consistent parameters. Note that Ps is 
found to be nearly T-independent below the pseudogap 
onset. Equation 16 is used to fix the value of u. From 
Eq. C44, Ps oc u~^, while the cr-model calculations^^'^^ 
give Ps = JS^. Equating these two expressions for a; = 0, 
T = gives = 0.256eV, which is assumed for all 
dopings. The calculated values of ps are illustrated in 
Fig. 15c, based on Eqs. C45, C44. 
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IV. ARPES SPECTRA 
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FIG. 15. (a) 5 — So (thin solid lines — —So); (b) ps calcu- 
lated from Eqs. C44, C38; (c) -So; (d) A, Eq. 22. In all the 
plots, the solid curves correspond to a; = 0.0, dotted lines: 
X = —0.04, short dashed lines: x = —0.10, long dashed lines: 
X = -0.15. 



A. SCR Transition and Correlation Length 



B. General Results 



Given the above parameters, the doping dependence 
of the MF and SCR transitions is compared in Fig. 15 
for the four electron dopings studied in Refs. 10, 11. 
The MF transition occurs when the bare Stoner factor 
Sf) = 1 — XqqU becomes negative. Fig. 15a. However, 
in SCR the renormalized Stoner factor S stays positive, 
so there is no T > phase transition (Mermin- Wagner 
theorem), although S — Sq has a strong increase near the 
temperature where Sq changes sign. There is still a zero- 
T Neel transition, controlled by the quantum corrected 
Stoner factor, Sq = ij — XqqU . From Fig. 15c, it can be 
seen that at x = —0.15, the system is close to a QCP, 
6o{T 0) ^ 0. This QCP is controlled by the Stoner 
criterion of the zero-T antiferromagnet. While there is 
no long range order, there is still a Mott (pseudo)gap, 
controlled by short-range order. Fig. 15d. This will be 
discussed further below. 

In the renormalized classical regime, the vanishing of 
(5 as T ^ is controlled by a correlation length, which 
can be written as^^ 



2-nps/kBT 



(16) 



with spin stiffness ps give n by Eq. C44. The prefactor 
^0 is T-dependent, = ^Ae/2CT. This T-dependence 
agrees with the one-loop cr-model results^^ rather than 
the more accurate two- loop results^*'®". This difference 



Given the susceptibility, Eq. 3, the self energy can be 
calculated as 



Y.{k,iujn) 



Go[k + q,iuJn + ii^m)Do{q,'i 



9^X0 
V 



E 



Ct 



with bare Green's function Go{k^iujn) = \/{i(jJn — S,^) 
and magnetic propagator Dq, Eq. C26; for the form of 
the integral, see the discussion near Eq. C33. In addition, 
Xo — XoiQi 0), (f = Q + (f , n is the Bose function, and 



g'^Xo = U^Xo{Uxo{Q,ii^n) + 



1 ) _ 3?/ 



(18) 

(Ref. 61). The last form is an approximation based on 
the empirical substitution xq 1/U in the pseudogap 
regime. After analytical continuation, the imaginary part 
of the retarded self energy is 



V 



I de[nie) + fiij:^,)]x 



x,5(c^ + e-^g ) 



Ce 



(<5 ^g'2)2 -h (Ce)2 ■ 
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The resulting self energy is plotted in Fig. 16 for T = 
IQQK . (The weak oscillations seen in some branches of 
S/ are an artifact due to an insufficient density of points 
in the numerical integration.) Note that /mE has the 
form of a broadened 5- function peaked at w = Cj+q- If 
it were a (5-function, /mE = —i:K^5{uj — Cj+g)? then 



1 



de 



A2 



(20) 



which is equivalent to Eq. 22. (In the last equality in 
Eq. 23 the limit — > is an adiabatic approximation^^, 
while the approximation is made that x peaks near Q.) 
Thus, as long as there is short-range magnetic order (A. 
or ps non-zero), there will he a Mott (pseudo)gap. Equa- 
tion 23 was also derived by Schmalian, et al.^^, but they 
did not discuss its significance. 



C. Application to the Cuprates 



SO away from the (5-function 
G(fc,cj) 



- - Re^^ik, uj) (^ - efe)(t^ - C^+q) - A2 

(21) 

This is exactly the Green's function of the mean field 
calculation^ with the substitution A — > A, where A 
can be evaluated by integrating 



A^ = -i 

TT 



dLoIm^^{k,uj) 



(22) 



Fig. 15d. This result is due to the Bose term n(e) in the 
square bracket of Eq. 19, the Fermi function / making 
no contribution. This leads to A being independent of k. 



200 
180 
160 
140 
120 
uj- 100 
30 
60 
40 
20 





(1,1) 








(1.0) 






(0,0.2) 






(0.5,0.5) 








(0.0.6) (0,0.4) 




\ 


I , X 







-3-2-10123466 
(j/t 



FIG. 16. Imaginary part of the self energy, Eq. 19, assum- 
ing 1/C = O.OSt, 5 = 0.002, a„ = 1, T = lOO.K'. The branches 
are labelled ik^ 



, ky), in units of tt. 



Equations 21,22 constitute an important result: the 
connection between the Mott gap and short-range mag- 
netic order. Recalhng that A = U < Mi >, or A^ = 
C/2 < S, >2, where < M, >= (-1)* < 5^ > is the 
staggered magnetization, then, in the spirit of an alloy 
analogy, a short-range order parameter can be defined as 



-g 

4/3 



J2 < S^+iT)SJ-iO) > e'^^dr 
<»,i> 

c,)x+_(fc,zu;)^|^^X+-(fc,0) (23) 

k 
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FIG. 17. Spectral functions for (a) x = 0, (b) a: = -0.04, 
(c) X = -0.10, and (d) x = -0.15, at T = WOK. Solid lines 
at (7r,0), and long dashed lines at {n/2,n/2). 

Using the correct ImY,^ from Eq. 19, and the calcu- 
lated parameter values from Table I, ARPES spectra are 
calculated for electron-doped cuprates, at the four dop- 
ings for which detailed data are available^". Figure 17 
shows typical calculated spectra for several fc-points in 
the a-b plane. Broadened Hubbard bands are found, 
which gradually smear out at high temperatures as S 
increases (^ decreases). There is a well defined pseu- 
dogap, with two peaks in the spectral function at a given 
k. It should be stressed that since there is no inter- 
layer coupling, long range antiferromagnetic order exists 
only at T = OK. The resulting dispersions are shown in 
Fig. 18. Figures 19- 21 illustrate the temperature depen- 
dence of Im{G) and /to(E), for two dopings, a; = and 
-0.15. The broadening of the peaks can be understood 
from Eq. 19: particle-hole excitations are present within 
a range -iza^/C of ^jr_|_^- Away from this particle-hole 
continuum the main peaks are sharp, while they broaden 
when they enter the continuum. 
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FIG. 18. Dispersion relations for electron doped materials, 
calculated at T = lOOK: (a) a; = {U/t = 6), (b) a; = -0.04 
(U/t = 5), (c) X = -0.10 {U/t = 3.5), and (d) x = -0.15 
{U/t — 2.9). Weaker features are denoted by smaller circles; 
for X — —0.15 all shadow features are extremely weak. 



20 
.-.15 

5 




(a) Si 



A! 




0.10 0.15 0.20 0.25 0.30 

(0 (cV) 




-0.5 -0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 
(i)(cV) 



FIG. 19. Temperature dependence of (a) spectral func- 
tion and (b) imaginary part of self energy, for x = 0.0 at 
(7r,0). Temperatures are 100, 500, 1000, 2000, 3000, 4000, 
and 5000K. 



FIG. 20. Temperature dependence of (a) spectral func- 
tion and (b) imaginary part of self energy, for x = 0.0 at 
(7r/2,7r/2). Temperatures are 100, 500, 1000, 2000, 3000, 
4000, and 5000K. 

In contrast, for x = —0.15, Fig. 20, the splittings are 
absent near {tt/2, n/2), and vanish near (tt, 0) by ~500K, 
and the lines actually sharpen on warming. If the effec- 
tive U is reduced to 2.5t, no splitting is found, but the 
peak position and broadening have an anomalous T de- 
pendence. Clearly, the system is very close to a QCP. 
Figure 22 shows in more detail how the spectrum evolves 
with U near this point. 



Note that the Mott gap collapse is anisotropic: for 
the undoped case, the nodal gap collapses between 2- 
3000K, while a gap persists near (7r,0) above 5000K. 
/m(S) has striking oscillatory structure, particularly 
near (7r/2,7r/2), which produces a similar weak struc- 
ture in Im{G) at low T. [Similar, weaker oscillations are 
present near (tt, 0), which can be better seen in Fig. 4c 
of Ref. 35.] In addition, there is a very intense, strongly 
T-dependent peak in Im{Yj) exactly at ^^r^^ (Fig- 19b - 
also present but not shown in Fig. 20b). It is the di- 
vergence of this peak as T ^ which signals the AFM 
transition. At low temperatures, the peak positions in 
Im{G) have a temperature dependence consistent with 
the collapse of the Mott gap - e.g., the LHB shifts to 
higher energies (toward midgap) at higher temperatures. 
Some experiments on hole doped cuprates find the oppo- 
site dependence^"^, which can possibly be understood as 
a localization or phase separation effect. 
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FIG. 21. Temperature dependence of spectral function for 
x = -0.15 at (7r,0), for U/t = 2.9 (a) and 2.5 (c). Tem- 
peratures are 100 (solid line), 500 (long-dashed line), 1000 
(short-dashed line), and 2000K dot-dashed line), (b): imagi- 
nary part of self energy at T = lOOK, U/t = 2.9. 
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FIG. 22. [/-dependence of spectral functions for x = —0.15 
at T = lOO-ft: near the T = QCP, for U/t = 2.5 (short dashed 
line), 2.7 (long dashed line), 2.9 (solid line), 3.0 (dot-dashed 
line), and 3.2 (dotted line). 
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FIG. 23. Temperature dependence of gap A for (from high- 
est to lowest) X = 0, -0.04, -0.10, and -0.15. Arrows show 
mean field transition temperature Tjv. 

Thus, the SCR calculation agrees with the mean-field 
results^^, if the mean-field gaps and transition tempera- 
tures are interpreted as the opening of a pseudogap at 
finite T, with the long-range AFM appearing only at 
T=0. A direct comparison of the transition temperatures 
is presented in Fig. 15d, presented on a linear T scale in 
Fig. 23. Moreover, the overall dispersions, Fig. 18 are in 
quite good agreement with the mean field results^'^ and 
experiments^*^. This agreement is somewhat surprising, 
since the model is not fully self-consistent. For instance, 
the parameter C involves Landau damping of the spin 
waves by electron-hole excitations, and hence depends on 
the electronic dispersion near the Fermi level. Thus, the 
opening of the pseudogap should have a strong influence 
on C, which is not accounted for. 

Finally, Fig. 24 displays Fermi surface map for x = 
—0.04 and —0.10. The pseudogap along the zone diag- 
onal associated with the hot-spot scattering^^ is clearly 
seen. These should be compared with the mean-field^^ 
and experimental^^ results. One interesting aside: in the 
mean field calculation, with sharply defined bands, it was 
necessary to include a t" parameter to reproduce the ex- 
perimental hole pocket near the zone diagonal. In the 
SCR calculation the spectral function peaks are consid- 
erably broader, and no t" parameter is needed. 
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FIG. 24. Fermi surface map for x = —0.04 (a) and —0.10 



(b). 



V. THREE DIMENSIONAL NEEL ORDER 

In the physical cuprates, the interlayer hopping has an 
anomalous dispersion, generally written as = tzoicx — 
CyY' . This formula holds for bilayer splitting, and in gen- 
eral when the Cu02 planes are stacked uniformly. How- 
ever, as explained in Appendix E, many of the cuprates, 
including NCCO, have a staggered layering.^ with the Cu 
in one Cu02 plane laying above a vacancy in the neigh- 
boring Cu02 sheet. This leads to a magnetic frustration: 
the Cu in one sheet has four nearest neighbors in the ad- 
jacent sheet, two with spin up, two with spin down. This 
frustration is reflected in a more complicated dispersion 
of t^: 



tz = t^icx - Cyf COS cos ^ 



a kya 
cos 
2 2 



(24) 



which vanishes at (tt, 0) and (0, tt), and leads to a greatly 
reduced interlayer coupling. (Effects of AFM frustration 
associated with layering have been discussed in Ref. 65.) 

The consequences of both uniform and staggered stack- 
ing are explored in Appendix E. If the c-axis resistivity 
is coherent, it can be used to estimate the interlayer hop- 
ping tzo- It is found that the value of t^Q needed to 
produce a given resistivity anisotropy is approximately 
5 times smaller for uniform stacking, to account for the 
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frustration in the staggered stacking. With the corre- 
sponding izo's determined from resistivity, both forms 
of interlayer coupHng give rise to comparable interlayer 
couphng, and hence a finite Neel temperature. While the 
optimal Q-vector depends on doping, at half filling both 
forms predict Q = (tt, tt, 0), consistent with experiment 
in La2Cu04. Even for quite strong anisotropy, this mech- 
anism can account for the observed T^rs (in fact, tends to 
overestimate Tat), without the necessity of invoking ad- 
ditional mechanisms, such as a Kosterlitz-Thouless tran- 
sition, with the reduced spin dimensionality caused by 
spin-orbit coupling effects^^"^^ . 

Within mode coupling theory™ (Appendix E), the 
Neel temperature is found from the gap equation 
(Eqs. El, C51) 



Xo{T)U = 7^ + 



SuTa^ In(^) 



ttA 



(25) 



where T^o ^ tl is defined below Eq. E8. It is found that 
T3D is approximately constant, independent of doping in 
the electron-doped regime. Apart from a small numerical 
factor, Eq. 25 diff'crs from the isotropic three-dimensional 
result by the logarithmic factor, which diverges (T/v 0) 
as tz — > 0. 

Equation 25 can be rewritten in a suggestive form. Ap- 
proximating ps by = A{x^U — ri)/12ua? (Eq. C44), 
then, using Eq. 16, the Neel transition occurs when 



N 



rr, 



N 



(26) 



where = J(^^o/^)^ J = 4:t^/U, and T = UIo/UTsd- 
A very similar form was proposed earlier^^. 




stacking, with adjusted to reproduce the observed re- 
sistivity anisotropy, but Appendix E shows that the over- 
estimate is generic: the coefficient of the logarithm needs 
to be larger to reduce T^. Also shown in Fig. 25 (dotted 
line) is a simplified model, which assumes that 



nA 



Sua^ ln( j?^) 



(27) 



is doping independent, Tq* — 1200iir. This model repro- 
duces qualitatively the shape of the numerical calcula- 
tion, but with a magnitude comparable to experiment. 
The magnitude of T/v could be matched almost quanti- 
tatively if Uef f also has a significant temperature depen- 
dence, as discussed in Appendix E . The overall doping 
dependence is also comparable to experiment. The agree- 
ment could be further improved by using a smaller value 
of t' , which would shrink the doping range over which 
Neel order occurs. 



VI. DISCUSSION 

A. Magnon Bose Condensation 

Figure 26 shows the sharp peak which arises in /mE at 
low T. The growth is exponential, approximately match- 
ing that of the coherence length, Eq. 16. (Note that it 
requires a fine mesh in the integral of Eq. 19 to capture 
this growth.) This peak arises exactly at the incipient 
magnetic zone boundary, and turns into true Bragg scat- 
tering at the transition to long range order: the increase 
in peak height is almost exactly compensated by a de- 
crease in the width of the peak. A simple physical expla- 
nation is that the SDW transition can be interpreted as a 
Bose condensation of the zone boundary magnons. Then 
the Mermin- Wagner theorem reduces to the fact that in 
a two-dimensional system, Bose particles can only con- 
dense at r = 0. A similar explanation for the transition 
has been presented earlier'^. 
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FIG. 25. Comparison of experimental Neel temperatures 
for NCCO and LSCO, (solid line), and for the stripe (mag- 
netic) ordering transitions observed in Nd-substituted LSCO 
[ 52] (solid line with squares) with the model of interlayer 
coupling with staggered stacking and tzo = t/10 ~ 30meV , 
plotted as Tjv/lO (dot-dot-dash line). Also included is the ap- 
proximate expression, Eq. 27 (dotted line with circles). (Note 
that there is a range of hole doping for which A is found to be 
negative; in this range Tjv was arbitrarily assumed to vanish 
in the staggered model, Tjv = 0.) 

Figure 25 compares the calculated value of Tjv with 
the experimental values. While the overall doping depen- 
dence is comparable, the calculated T/v is about an order 
of magnitude higher. The calculation is for staggered 
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FIG. 26. (a,b): Blowups of 7m(E) for x = at (tt, 0) (a) 
and (7r/2,7r/2) (b) at T = lOOK (solid lines), 500K (long 
dashed lines), and lOOOK (a) or 750K (b) (short dashed 
lines), (c) Maximum of Jm(E) vs T for (tt, 0) (squares) 
and (7r/2,7r/2) (circles); 0.1/(full width at half maximum) 
for (7r,0) (triangles) and (7r/2,7r/2) (diamonds); solid line = 
corresponding £,{T), Eq. 16. 

B. SCR Calculation of NAFL Parameters 

The SCR is perhaps the simplest model in which fluc- 
tuation effects are included to satisfy the Mermin- Wagner 
theorem, allowing one to ask questions such as, how does 
a Mott gap appear if there is no long-range order to gen- 
erate a smaller Brillouin zone? It can be seen that the 
Mott gap is really a pseudogap, and is relatively insen- 
sitive to the appearence of long-range Neel order. How- 
ever, it must be kept in mind that the original SCR^^ 
is intended to describe weak itinerant ferromagnets and 
AFM's, and will have to be extended to account for the 
strong renormalization of the electronic bands. 

Given these limitations, the present paper attempts 
to use the SCR to calculate the parameters. A, C, etc., 
for an NAFL-like model of the low energy physics. The 
principal findings include: 

• the program fails for hole doped cuprates, since the 
parameter A is negative, signalling an instability of the 
commensurate AFM phase against either incommensu- 
rate SDW phases or nanoscale phase separation. 

• For electron doped cuprates, considerable progress 
can be made, but the model still has problems, in partic- 
ular a poorly defined spin-wave interaction parameter u. 



Here u is estimated by comparison with sigma model re- 
sults near half filling, but its possible doping dependence 
is unknown. 

• The present theory differs from conventional NAFL 
theory by the inclusion of two cutoff parameters, qc and 
w^. These cutoffs shrink to zero at either the H- or C- 
points, and in particular cause the A parameter to have 
a strong temperature dependence in the electron-doping 
regime. 

• The combined SCR and mean-field^^ results are in 
excellent agreement with the ARPES data^*^ on the dop- 
ing dependence of the Mott transition in the electron- 
doped cuprates. 

• While the collapse of the Mott gap and the termina- 
tion of the hot spot regime are in principle independent 
QCPs, in practice they fall very close to each other. The 
present calculations and Ref. 20 explain why this is so. 

These points are discussed in more detail in the follow- 
ing subsections. 

C. Mott Transition vs Slater Antiferromagnetism 

In reviewing the history of magnetism, Anderson re- 
cently pointed out^^ that the theories fall into two dia- 
metrically opposed classes: band theory and atomic mod- 
els, and the latter are typically more successful. In partic- 
ular, for strongly correlated Mott insulators Mott's model 
of nearly localized spins seems like a better starting point 
than Slater's theory of spin density wave antiferromag- 
netism. The main problem is that the Slater theory ties 
the Mott gap to antiferromagnetic order, predicting too 
high a Neel temperature. 

While a localized picture may be a more convenient 
starting point near half filling, nevertheless it should be 
possible to develop a (perhaps more complicated) picture 
based on band theory. There is a general desideratum to 
be able to extend band structure calculations to all mate- 
rials; moreover, this is important in the present instance 
because the doping dependence of U suggests that there 
is a crossover from strong to intermediate correlations, 
which may be better handled by working throughout in 
a band structure formalism. The present mode coupling 
calculation seems to be an appropriate starting point. 
The Mermin- Wagner theorem leads to a decoupling of the 
Mott gap and Neel order. How this might be extended 
to three-dimensional Mott insulators remains unclear. 

Local physics should show up in a band structure 
calculation as a very narrow band width. In a mean 
field SDW calculation, the width of the Hubbard bands, 
^ t^/U ^ J, is not small. 
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D. Stoner Criterion and Crossover from Small to 
Large FS 

In the Hartrcc-Fock model, the crossover from small to 
large Fermi surfaee is coincident with Mott gap eoUapse, 
and comparison with experiment suggests that this cor- 
rectly describes the situation in electron doped NCCO. 
The SCR theory confirms this result and offers additional 
insights. Fluctuations preclude long-range order at finite 
temperatures (in the absence of interlayer coupling) , so 
the physics is that of a zero temperature AFM QCP, con- 
trolled by a Stoner criterion: the QCP occurs when x/7 
is too small. At finite temperatures, the Mott gap is 
replaced by a pseudogap, and the FS crossover still co- 
incides with the collapse of the Mott pseudogap, which 
occurs at the QCP. Whether a similar crossover can oc- 
cur in the absence of Mott gap collapse remains to be 
seen. 



E. Hole Doped Cuprates 

The present calculations make two predictions for hole- 
doped cuprates: (1) there are strong indications for insta- 
bility against phase separation and stripe physics; but (2) 
despite this, the termination of strong magnetic fluctua- 
tions should be approximately electron-hole symmetric^'' 
- associated with the symmetric susceptibility plateau in 
Fig. la. 

1. Stripes 

The situation in hole doped cuprates is complicated by 

the presence of stripes. At all levels, mean-field, RPA, 
SCR, striking differences between hole and electron dop- 
ing are clearly revealed. All techniques provide strong 
evidence for the instability of the AFM state for hole 
doping, while it remains stable under electron doping. 
In earlier Hartree-Fock and RPA calculations, evidence 
for the instability of the hole doped phase was found in 
the ordered magnetic phase, within the smaller magnetic 
Brillouin zone: e.g., the spin wave dispersion is unsta- 
ble. In the SCR approach, there is no phase transition at 
finite temperatures, but even in the paramagnetic phase 
there is evidence for the instability. A detailed analysis of 
the real part of the bare susceptibility (as appropriate for 
a Stoner criterion) provides evidence for instability of a 
commensurate magnetic phase at Q (negative curvature 
A). 

2. Pseudogap 

In hole doped cuprates, ARPES finds two features 
which are commonly referred to as pseudogaps - a 'hump' 
feature found near (tt, 0) at higher binding energy than 



the main, superconducting 'peak', and the 'leading edge 
gap', a loss of spectral weight in the immediate vicinity 
of the Fermi level. This latter feature is not explained 
by the present calculation; it may be associated with the 
onset of strong superconducting fluctuations^^''^^. Alter- 
natively, such a gap has been found in a dynamical cluster 
expansion calculation of the Hubbard model 

On the other hand, the 'hump' feature can be 
consistently interpreted as the collapse of the Mott 
pseudogap^^'^". While ARPES only sees the feature be- 
low the Fermi level associated with the lower Hubbard 
band (LHB), tunneling ''"^'^''^'''^ finds an approximately 
symmetrical peak feature above the Fermi level, asso- 
ciated with the UHB - as if the pseudogap were pinned 
to the Fermi level. This has led to a number of alter- 
native models for the pseudogap, in terms of supercon- 
ducting or charge density wave (CDW) fluctuations (the 
latter possibly related to stripe physics). Recently^" it 
was found that the same mean field model of Mott gap 
collapse can approximately explain the data (see Fig. 2b). 
In this model, the lower pseudogap peak is the VHS of 
the LHB, while the upper peak is due to the leading edge 
of the UHB. As the Mott gap collapses, the two features 
merge. A careful tunneling study of the 'hump' features 
could look for the predicted asymmetry of the features 
about the Fermi level. Such a study should best be done 
in single layer Bi2Sr2Cu06,^^ where complications due to 
bilayer splitting are absent. 

The above interpretation requires that for hole doping 
also the Mott gap must collapse slightly above optimal 
doping. This is consistent with recent experimental ob- 
servations of a QCP21. The early SCR results of Table II 
point in the same direction: the very small and positive 
values found for the Stoner parameter yo in optimally and 
overdoped cuprates suggests the proximity to a QCP near 
optimal doping. Moreover, the model predicts that at the 
QCP, where the pseudogap just closes, the Fermi level is 
exactly at the VHS (H-point). This result had been found 
experimentally in some lightly overdoped cuprates^^'^^. 

While the doping of Mott gap collapse is approximately 
electron-hole symmetric, some significant differences re- 
main. Thus, for electron doping the magnetic correlation 
length remains large up to the QCP, while the correla- 
tion length is only a few lattice constants on the hole 
doped side. This may be related to competing order - 
indeed, in Nd substituted LSCO, long-range incommen- 
surate magnetic order is found^ up to a; = 0.2. Alterna- 
tively, Schmalian, et al.^^ were able to reproduce a hump- 
like pseudogap with small correlation length, by a careful 
summation of the full diagrammatic perturbation series. 
[Note that near a VHS, all competing electron- hole insta- 
bilities - SDW, CDW, flux phase, shear ('Pomeranchuk') 
instability will lead to similar pseudogaps near (7r,0), 
and indeed in the presence of strong fluctuations, all will 
contribute in a comparable fashion, A*^ Af .] 
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F. VHS 

1. Electron-Hole Asymmetry 

To study the role of the VHS in the Mott transition 
and high- Tc superconductivity, one would ideally like to 
study a system in which one could turn the VHS on or 
off. Electron vs hole doping of the cuprates would appear 
to approach this ideal. In switching from electron to 
hole doping (from NCCO to LSCO) Tc increases by less 
than a factor of two, and apparently remains d-wave, 
while the normal state properties change drastically, with 
nanoscale phase separation on the hole doped side only. 

On the other hand, the one hole vs one electron sys- 
tems should be much more similar: in either case, nearest 
neighbor hopping is frustrated by breaking local antifer- 
romagnetic order. Hence in both cases, the low energy 
states will be magnetic polarons. Indeed, the electron- 
doped polarons may be more localized, since there is no 
interpolaron attraction (i.e., tendency to phase separa- 
tion). The ARPES spectra for low-electron doping {x=- 
0.04) show an additional pseudogap at the Fermi level, 
which may be related to localization. 

2. Temperature Dependent VHS 

As noted by Onufrieva and Pfeuty''", the VHSs associ- 
ated with the susceptibilities (and hence with charge or 
spin nesting) are different from those associated with the 
density of states (and superconductivity). Thus, whereas 
superconductivity will occur at the same optimal doping 
for all temperatures, the doping of maximal nesting in- 
stability is a strong function of temperature. 

This contrasting behavior of nesting vs pairing suscep- 
tibilities is related to a characteristic difference in the 
nature of the two instabilities. A superconducting in- 
stability has an intrinsic electron-hole symmetry, which 
means that the gap is tied to the Fermi level, and a full 
(s- or d- wave) gap can be opened at any doping level. 
On the other hand, a nesting gap is dispersive, and only 
part of it lies at the Fermi level (except in special cases). 
Furthermore, a (superlattice) Luttinger's theorem must 
be obeyed, requiring the presence of residual Fermi sur- 
face pockets. Stated differently, a full nesting gap can 
only open at integer filling, so as the interaction strength 
increases, any nesting instability must migrate to inte- 
gral doping (e.g., half filling in the original band struc- 
ture). This same VHS migration is mirrored in the T- 
dependence of the magnetic (or charge) susceptibility. 

3. VHS Transitions 

We have seen that the doping-dependent f/g// gives 
rise to a Mott gap collapse near the edges of the suscep- 
tibility plateau in Fig. 1. If U^ff is smaller (dot-dashed 



line: Ueff reduced by 2/3), more complicated behavior 
should arise. Due to the peak in x near the H-point, 
there could be a reentrant transition, with one magnetic 
order near half filling, and a second near the VHS. For an 
even smaller U^ff (or replacing U^ff — » J)^^, the transi- 
tion near a; = can be eliminated, leaving a spin density 
wave transition near the VHS. In principle there could 
even be a phase separation between two AFM phases: an 
insulating phase near half filling and a metallic phase 
near the VHS. 



G. Future Directions 

It must be stressed, however, that the present theory 
is not fully self-contained. There are three significant 
limitations. First, an improved calculation of the doping 
dependence of the Hubbard U parameter is a desidera- 
tum, perhaps along the lines of earlier calculations^"'*^'^^. 
Since all of these calculations lead to different doping de- 
pendences for U{x), the actual doping dependence must 
be regarded as an unresolved issue. 

Second, several unsuccessful attempts have been made 
to calculate the quartic interaction parameter u. A sim- 
ple ansatz for u{x) is introduced, based on consistency 
with the t — J model, which leads to good results. A 
deeper understanding of why this works, and whether u 
is doping dependent, is desirable. 

Finally, the susceptibility was calculated with the bare 
electronic bands, but when the bands are renormalized 
to first order in x a pseudogap opens. This gap should 
be self consistently incorporated into the calculation of 
X, as in the FLEX and spin fermion approaches; it is 
expected to have a profound effect on the temperature 
dependence of the parameters (especially A) and on the 
residual density of states in the pseudogap, particularly 
near a; = 0. 



VII. CONCLUSIONS 

The main results of this rather long paper are briefly 

summarized: 

• Fluctuation effects were added to the mean field Hub- 
bard model via a mode coupling calculation, which al- 
lowed satisfying of the Mermin- Wagner theorem (Tjv = 
0). It was found that the mean- field gap A^/ and Neel 
temperature T'^^ evolved into a pseudogap Ap^ A^/ 
and an onset temperature T* ^ T™^ (as is familiar from 
the related CDW results). 

• The resulting dispersions and Fermi surfaces are in 
excellent agreement with photoemission experiments on 
electron-doped cuprates^", while the pseudogap seems 
consistent with ARPES and tunneling results in hole 
doped cuprates^°. It is interesting to note that a re- 
cent t — t' — t" — J model calculation seems consistent 
with the first doped carriers forming weakly interacting 
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quasiparticles in pockets of the respective upper or lower 
Hubbard bands, for either electron or hole doping®^. 

• The zero-temperature Neel transition is controlled by 
a Stoner-like criterion, hence is sensitive to the bare sus- 
ceptibility and in turn to the Fermi surface geometry (hot 
spots). This lead to an approximately electron- hole sym- 
metric QCP near optimal doping (termination of hot spot 
regime), at which both zero temperature Neel transition 
and pseudogap transition simultaneously terminate. 

• The model leads to a NAFL-type susceptibility, and 
the calculation of the NAFL parameters has been re- 
duced to a calculation of the coupling parameters U and 
u, the former having a significant doping (and possibly 
temperature) dependence. At present, U{x) is estimated 
from experiment, and the mode coupling u via consis- 
tency with the t — J model. (A small portion of the 
renormalization of U arises from quantum corrections to 
the Stoner criterion.) 

• Whereas the antiferromagnetic state at Q is stable to 
electron doping, hole doping leads to an incommensura- 
bility, which is interpreted as an indication of instability 
to phase separation (as found in the mean field calcula- 
tions). This asymmetry follows from the properties of 
the VHS. 

• Finally, a striking temperature/frequency dependence of 
the susceptibility peak, from the VHS at low T to half 
filling at high T, found earlier^'', is interpreted in terms 
of Luttinger's theorem: if the coupling is strong enough 
to open a full gap, the gap must fall at half filling. 
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APPENDIX A: THREE BAND MODEL 

A major simplification of the present calculation is to 
treat the cuprates in a one-band model. This is con- 
sistent with the Zhang- Rice picture^^, although the ap- 
proximation is less drastic for electron doping, since the 
upper Hubbard band is already predominantly copper- 
like. Nevertheless, the model also describes the doping 
dependence of the 'lower Hubbard band', which is really 
a charge transfer, predominantly oxygen-like band. Here 
an explanation for why this simplification works is sug- 
gested. 

Even without carrying out self-consistent calculations, 
the nature of the Mott transition can be understood by 
introducing a doping dependent gap. The energy bands 
can be calculated from the hamiltonian matrix 



H = Y^ Ad]dj + tcuo[d]p^ + (c.c.)] 

j <i,j> 

+ too[p]pr + (cc.)] (Al) 

<jj'> 

where A is the difference in on-site energy between cop- 
per and oxygen, tcuO is the copper-oxygen hopping pa- 
rameter, and too the oxygen-oxygen hopping parameter. 
For good agreement with the doping dependence of the 
one band model, it is necessary to properly incorporate 
the Hartree correction to the self energy, A = Aq + T,h, 
J^H — Uui (for up spins), and = n/2 — mg, with n the 
average electron energy. In Figs. 27, Aq = is assumed. 

The band dispersion is extremely similar to that found 
in the one band model. Fig. 3 of Ref. 11, even though 
the lower band crosses over from the Zhang-Rice (hy- 
bridized copper-oxygen band) at half filling to a more 
copper like lower Hubbard band with increasing electron 
doping. In addition, the effective magnetizations are pro- 
portional (inset. Fig. 27d), although the one-band model 
overestimates the magnetization by 1/3. This can be 
understood: in the three-band model, the shape of the 
Hubbard bands is fixed by the combined effects of the 
magnetic instability and hybridization with the oxygen 
band. In the one band model, only the former effect is 
present, necessitating a larger value of m to produce the 
same net splitting. 

This remarkable agreement between one and three- 
band models goes well beyond the Zhang-Rice model. 
That model is restricted to the LHB in a small range of 
doping near half filling; the present results compare both 
LHB and UHB over the full range of electron doping. 
The result is nontrivial - in the three band model, the 
bonding and non-bonding bands are also split into up- 
per and lower Hubbard bands. This degree of agreement 
comes about because the parameter A is approximated 
by the two components (magnetic and nonmagnetic) of 
the Hartree term. In turn, this suggests that in the ab- 
sence of magnetic effects the Cu and O energies are nearly 
degenerate - as found in early LDA band structure cal- 
culations (see discussion in Ref. 84). 
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FIG. 27. Dispersion of two antibonding bands in 
three-band model, assuming mq = 0.3 (a), 0.2 (b), 0.05 
(c), and 0.01 (d). Inset to d: effective magnetization 
rrieff — mU/Qt for the three-band (solid line) and one-band 
(filled circles) models. The one-band result has been multi- 
plied by 3/4 to better agree with the three-band results. 



APPENDIX B: CHARGE SUSCEPTIBILITY AND 

Ueff 

The proper choice of vertex corrections is an unre- 
solved issue in the analysis of the Hubbard model. It 
is known to be of critical importance for generating a 
pseudogap^^. Here, by comparing simple mean- field and 
SCR models to experiment, it is shown that the net ef- 
fect of vertex corrections is to make the coupling U ef- 
fectively doping (and possibly temperature) dependent. 
Kanamori"*^ showed that the effective Hubbard U should 
decrease with doping, as an electron can hop around, 
and hence avoid, a second electron. In the limit of a 
nearly empty (or full) band, this should lead to a correc- 
tion of the form U^ff ^UI{1 + U/W), where W = U 
is the bandwidth. It was found®^'^'^ that Monte Carlo 
calculations of the susceptibility of a doped Mott insula- 
tor were approximately equal to the RPA susceptibility 
with suitable t/e//, and Chen, et al.^^ suggested the ex- 
plicit form Ueff = U/{1+ < P > U), with P given by 
a vertex correction to the susceptibility and < • • • > an 
average over q, at zero frequency. Figure 28b presents 
a calculation for Ueff based on Chen, et al. However, 
whereas Chen, et al. performed the average in the para- 
magnetic phase, using bare Green's functions, here the 
dressed Green's functions appropriate to the Neel phase 
are used, to approximately incorporate the effect of this 
gap. This makes little difference, since P is dominated 
by the intraband terms, and hence remains finite at half 
filling. Explicitly, 



P = -1 ^ J7,,, (fc, k + q)hj (fc, fc + q). 



E^{k) + Ej{k') - w - iS 



E±{k) = -(e/c + Efc+g ± Eq), 



Ea — J (efc — €k+q ) 



-4A2 



(Bl) 



(B2) 



(B3) 



(B4) 



U^,j{k,k') = - {l + iAk){l+jAk') + ijBkBk', (B5) 

with i,j summed over — , A the AFM gap, and Ak — 
{ek~ek+Q)/Eok, Bk = A/Eok- In agreement with Chen, 



et al., the calculation finds U to be renormalized by a 
factor of 2 at finite doping, but does not recover a large 
U near half filling, although different results are found 
depending on whether a; = from the start (triangle) or 
whether x ~* from the hole or electron doping sides. 
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FIG. 28. Calculated Ueff assuming (a) simple screening or 
(b) full vertex correction of Chen, et al. [ 81]. In both cases, 
a bare U — 6.75t was assumed. Solid lines — electron doping; 
long dashed lines = hole doping; triangles (squares) in a = 
paramagnetic screening of U, at T = IK (2000K); triangle in 
b = undoped; circles = data of Ref. [ 10]. 

For modelling purposes, it is useful to have a Ueff 
which evolves smoothly from a large value at half filling to 
a reduced, Kanemori value at finite doping. A simple toy 
model consists of taking the RPA screening of a charge 
response. There should be a close connection between the 
Kanemori mechanism and screening. Screening involves 
creation of a correlation hole about a given charge, while 
Kanemori's Uef / involves the ability of a second charge to 
move around the first, while avoiding double occupancy. 
Near half filling, the second charge must move in the 
correlation hole. Approximating^^ the vertex correction 
by the RPA screening of the charge susceptibility. 



u 



1+ < X > [/' 



(B6) 



it is possible to reproduce^^ the experimentally 
observed^" doping dependence, while matching the cal- 
culation of Chen, et al. away from half filling. Fig. 28a. 

In this calculation, issues of self-consistency are also 
important. To minimize screening at half filling, it is nec- 
essary to reproduce the gap in the susceptibility. Hence, 
the susceptibility in Eq. B6 is approximated by the charge 
susceptibility in the AFM state, Xo" from Eq. 2.24 of 
Ref. 1, evaluated with the bare U — 6.75t. (In princi- 
ple, at finite doping there is a coupling to the longitu- 
dinal magnetic susceptibility^^, but this is neglected for 
simplicity.) The importance of using the AFM suscep- 
tibility is illustrated in Fig. 28a: the solid and dashed 
lines show Ueff calculated using the charge susceptibil- 
ity in the Neel state, while the corresponding lines with 
triangles use the paramagnetic susceptibility at low T. 
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The latter calculation finds a nearly doping independent, 
but small Ueff] the former reproduces a large, weakly 
screened U near half filling. Such a difference is expected 
in terms of screening: when there is no gap at half fill- 
ing, the enhanced susceptibility should be better able to 
screen U, resulting in a smaller Uef / . This suggests that 
Ueff should have an important temperature dependence 
as the gap decreases - which in turn will cause the gap to 
close at a lower temperature. Figure 28a also shows that 
there is a weak temperature dependence of the screening. 
The calculations suggest that the large values of U found 
in the cuprates are characteristic mainly of the half filled 
regime and relatively low temperatures. A similar but 
larger screening effect was recently reported by Esirgen, 
et al.80. 

This procedure is still not fully self consistent. If there 
is a large difference between the bare U and the screened 
Ueff, the gap in x should depend on the actual Ueff. 
However, since Ueff ~ f7 at half filling, any simple im- 
provement will not significantly change the overall dop- 
ing dependence. This is the same kind of lack of self- 
consistency found for the SCR approach, and will be here 
neglected. 



APPENDIX C: PATH INTEGRAL 
CALCULATION 

1. Formalism 

The partition function of the Hubbard model can be 
written as a path integral*^: 

Z = J DC'fDCexp [- drL] , (CI) 
with Lagrangian 



/ 



Z = / £»0e- 



(C5) 



U u 

(C6) 

Fourier transforming in space and (imaginary) time, the 
trace term can be rewritten 

Tr ln{^Go' + V)^Tr ln{-G^^) - ^ -Tr{GoVY , 



(C7) 



with 



Go fc, icon) = — C8 



V = ^^=^(j)ik-k',iujn-i(^m.), (C9) 

with — In Eq.C7, the odd terms in n average to 

zero, so expanding the action to fourth order in ^ yields 

5 = i ^ Il2{q,iu!n}(p{q,iwn)4>{-Q,-i^n) 

^ 4:{l3Nof ^ 'n4(9i, iwi)<^(gi, za;i)^(^2, «W2) x 

X(/)((f3, 2^^3)0(94, icj4), (CIO) 

where the prime in the second sum means summing over 
all qi, u)i, such that qi = 0, = 0, 

li2{q,iL0n) = ^[l-Uxoiq,ii^n)], (Cll) 



(C2) 

The quartic term can be decoupled by a Hubbard- 
Stratonovich transformation 

Z = j D(f>DC^DCexp[- dTL{4>, C^,C)] , (C3) 

m ct, c) = ^ clidr - M)a. - E 

+ TE'^? + ^E'^'("'T-"a) (C4) 



(neglecting a term involving ni| + nti). Integrating out 
the Fermion fields leaves an effective action in the field 



Il4{q,ioJn) = X] Go{k,ien)Go{k + qi,ien + iioi) x 

xGo{k + qi+ q2,ien + iuii + iu>2)Go{k - q4,ien - iuJi), (C12) 
and xo is the dynamic susceptibility 

Xo(g,«w„) = Go{k,itri)Go{k + q,ien + ii^n) 

PiVo z — 
k.ien 



j_^.fik+i_mi 



(C13) 



The mean field solution corresponds to assuming 

(t>i = ^oe'^^-^i and finding the saddle point solution of 
Eq. C6. Including interactions by summing bubble or 
ladder diagrams^ leads to the RPA susceptibility (see 
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Eq. C14 below), from which the spin wave spectra are 
calculated^'^. Fhictuations about the moan field solution 
are described by Eq. CIO. Due to the Mermin- Wagner 
theorem, these fluctuations are in the high temperature 
limit T > Tn. A naive perturbational analysis diverges 
(as demonstrated below), so a self-consistent analysis is 
necessary. Following the self-consistent renormalization 
(SCR) model of Moriya^^'^^, the exact dynamical suscep- 
tibility can be written as 



(C14) 



with the RPA susceptibility given by Eq. C14 with Xu = 
0. 

Solving Eq. C14 requires an equation for Ac/. An ap- 
proximate solution is found by replacing Xu{q, ioJn) by a 
constant 



A = Ac/(0,0) = xo( 



d^AF{M, T) 



)m=o, (C15) 



where the total free energy is written as F{M, T) = 
Fhf{M, T)+AF{M, T), with Fhf the Hartree-Fock free 
energy and 

AF{M, T) = -T V / dU[x{q, icon) - Xo{q, i<^n)] 
Jo 

= T ^[Zn(l - Uxo{q, ioJn) + A) + Uxo{q, iw„)]. (C16) 

q.n 

Equation C16 can be solved by expanding about the 
expected ordered state. The ordered states are found 
from the zeroes of the denominator of the dynamical sus- 
ceptibility, Eq. C14. For the present case, the largest 
bare susceptibility corresponds to antiferromagnetic or- 
der, q = Q. Then, defining 



S=l-Uxo{Q,0) + X, 



(C17) 



and So = 5 — X, it will be shown that ^ > 0, and ^ — > 

as T ^ that is, there is no finite temperature phase 
transition (the Mermin- Wagner theorem is satisfied). 



2. SCR Analysis 

Following the conventional analysis, Eq. C16 is ex- 
panded in terms of the small parameters u and (f = q— Q 
(analytically continuing iiVn co + ie): 

1-Uxo{q,co)+X = d + Aq'^ - Bw^ _ ^^^^ (^18) 
where the expansion coefficients are 

k 



r y^k+Q ^k> 



y^k+Q ^k> 



^ (^k+Q ^k) 
1 7f' 

-n^k)iT-.+ 



+ 



^^k+Q ^k^ 



16 (eg+Q - ej) {e^^^ - e^)^ 
-f"{e^)it'{l-c,Cy) + ^ -— ^)], (C20) 

y'^k+Q ^k' 

B = (C21) 



^^k+Q ^fe)^ 



C = -2nUY,[nei)5{^UQ-'k)\- 



(C22) 



These are similar to results found previously^^, special- 
ized to the dispersion of Eq. 1. However, with this dis- 
persion the A and B integrals formally diverge. Hence a 
more careful analysis is needed, presented in Appendix 
D. For now, A, B, and C will be treated as parameters. 



a. Free Energy 

Here, A is determined by minimizing the free energy, 
including quartic, corrections^^. With the varia- 

tional estimate 



F = Fo + -<S-So >So, 
the action S (Eq. CIO) becomes 



(C23) 



^ = -J (<^o + Aq'"^ + C|a;„|)(/)(g,ia;„)0(-g, -iUn) 



+uj dT j d^r^0(r,T)]^ (C24) 

and 5*0 is given by the same equation, with 5^ replaced by 
the variational parameter 5. Here u = li^a'^ /A(}Nq, with 
d = 2 and H4 evaluated at = (0,Q, 0, — Q), = 0, 
since higher order corrections are irrelevant. Then to 
Gaussian order 



<S- So >So= J ^{So-S) < (j){q,iun)'t){-q, -iuJn) >So 

(C25) 



< 4>{q,iuj„)(l}{-q,-iuJn) >So = 
Zq^ j dRe(f>{q, iu!n)dlm(j){q, iojn) x 

x0(g,iw„)0(-g,-iw„)e-5^o"'(9'^'^")"^(9''"")"^(-«'-^'^") 

dlniZo) 



= Do{q,iuJn) (C26) 
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with b^\q,iun) = iU/2)Da\q,iLUn) = {U/2)[S+Aq'^ + 
C|w„|]. Similarly, to Gaussian level 

Fo = —InZo = ^Y1 ^"^0 '(9. i^n), (C27) 

up to a constant. Writing the quartic term as Sint, the 
full partition function is approximated by Z ~ Zo[l— < 

Sint >Gauss +5 < Sf^^ >Gauss\, with 

< Sint >Gauss= U j dr J (ff< 4>{r,TY >Gausa 

= 3u^ dT j (fr[<(j){f,Tf >f, (C28) 

and 

< (t^{r,Tf >Gauss= ^ < (t){q,ii^n)(p{-q,-i'^n) > So 



^ J2 Do{q,iuJn)- 



[The term in <C S*^^^ ^ Gauss 

will not be needed.! 
The free energy Eq. C23 can thus be rewritten as 



q.lLOn 



^ m (^0 - 5)Do{q,iuJn) 



2/3 i 

+^[E D,{q,i^n)f 



(C30) 



with u = ni/NopV^ 



b. Stoner Factor 



The variational parameter is found from dF/d6 = 0, 



or 



(C31) 



The next step is to carry out the sum over Matsubara 
frequencies and wave vectors. For the former, using 



1 1 



ImX{e + i6) 



r°° e e 
■ / d-coth — ImX{e + id), (C32) 

Jo 2r 



then 



fd^qa^ r-'"^ de , 

/ ,^ ^o / COth 



e Ce 



2T {S + Aq'Y + (Ce)2 



. (C33) 



Note the sharp energy cutoff in Eq. C33. This comes 
about because the lincar-in-u' dissipation is a result of 
Landau damping of the spin waves by electrons near the 
hot spots, and therefore the dissipation cuts off when 
the spin wave spectrum gets out of the electron-hole con- 
tinuum. Numerical calculations (Fig. 11) show that the 
cutoff can be quite sharp, particularly near the VHS. 

Equations C31, C33 can easily be solved in the limit 
T = 0. In this case, there is a transition at 



So = -12m 



/2 2 ra^lC 



Jo 47r Jo 



de 



Ce 



"iuq^a^ 



-Ro = l-r], (C34) 



Ro 



1mi. 



tan ^{ttq) 



(C35) 



,q — Aq^/a^^. Since the right-hand side is finite 
and negative, fluctuations reduce but in general do not 
eliminate the order at T = 0. The quantum corrected 
Stoner criterion is Uxn = f], where representative values 
of T] are listed in Table I. 

However, for finite T, there are corrections ^ ln{8), so 
5 cannot be set to zero, and there is no finite temperature 
transition (the Mermin- Wagner theorem is satisfied). To 
see this, it is adequate to approximate coth(x) as 1/x for 
X < 1 and 1 for a; > 1. In this case, Eq. C31 can be solved 
exactly, Appendix A3. However, this exact solution is not 
very illuminating, and a simpler approximate solution 
will be given here. Since only the term proportional to T 
is singular, T and 6 can be set to zero in the remaining 
term. Defining 



So = So + r] - 



(C36) 



Eq. C31 becomes 
S-5o 



3uTa^, ,2CT 



/ -^tan-'{ 

T^^A y ' y ' 



ttA 



-ln{- 



S 



(C37) 



where the second line uses Eq. C50, below. Hence, there 
is no finite temperature phase transition, and S only ap- 
proaches zero asymptotically as T — > 0: approximately. 



5 = 2CTe~'^^'''°'/^"'^"^ 



c. Susceptibility 



(C38) 



The resulting susceptibility has NAFL form, Eq. 3, 
with (Eqs. C14,C18) the following explicit expressions: 



- ^ 



(C39) 
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C 



(C40) 



(C41) 



(C42) 



If the correlation length ^ is written in the form Eq. 16, 
and Eq. C51 is numerically solved for 5, then ps is exactly 
given by 



ksT . .A 



(C43) 



with ^0 = \J Using Eq. C38, an approximate pa is: 



„a ^1^0 1 



12ua2 



(C44) 



Ps is plotted in Fig. 15b, with u ^ = 0.256e\^, chosen to 
give a Ps in agreement with the results of Chakravarty, 
et al.^^ for a; = 0, T = 0. 



Io{x) = [ ln{tane)de = L{x) + L( J - x) - L( J), 
Jo 2 2 

(C49) 

and L{x) = — J^^ ln{cost)dt is the Lobachevskiy 
function^®. 

For most purposes, it can be assumed that S «t « 
Aq'^, a^, in which case Io{tan~^{x)) = 9{ln. (9) — 1), with 
6 = min{x, 1/x}, and then F2, Eq. C47, simplifies. 



.t. 



F2 = ln{-)[5 + Uan-\-)]+5- — 

Defining Z=l + {^ua? / 'k'^ AC)ln{a^ / 1) , then 

7K X 3Ma^^ , .2CT 
ttA 
which agrees with Eq. C37 when Z ^ 1. 



(C50) 



(C51) 



APPENDIX D: PARAMETER EVALUATIONS 



1. Hot Spots 



3. 'Exact' Solution of Eq. C31 

Approximating coth{x) = max{l/x,l), and introduc- 
ing the notation Aq"^ = Aq^ + 5, Ciq = Aq^/a^, and 
t = 2TC, the solution becomes 



■ [Fi + F2] , 



with 



rS+Aq- 

/ / 

Js Jo 



F2=t 



dx 



x'^ +y'^ 



t / —tan~^i- 

Js y y 



with 



(C45) 



au,tan~^ (dq) - - aojtan~^{ — ) (C46) 



(C47) 



h{x) = Ioitan~^{x)) - tan~^{x)ln{x), (C48) 



Hot spots are defined as the intersection of the Fermi 

surface with the line + ky = ir/a. At these points 
there is strong scattering, since the vector Q connects 
two hot spots. The hot spots dominate the integrals A, 
B, Eqs. C19,C21 at T = 0. The Fermi functions limit the 
integral to a sum of approximately wedge-shaped areas 
centered on the hot spots. Fig. 12. In the main text, 
the susceptibilities were calculated numerically. Here, 
analytical approximations are introduced to clarify the 
role of the hot spots. 

On a single wedge, a typical integral for, e.g., A be- 
comes 



/ 



d^q 



1 



1 



d^q 



1 



Letting ki = kto — k[, i = x,y, then to lowest order the 
energy difference becomes 



= ask'{l + Pgk') 
aQk±{l + l3ok\\), 



(D2) 



with ao = ao{sm6 + cos 6), (3$ = /3o(sin0 — cos 9), 
ao = 4:taSxo, and /3o = acxo/^s^o', k\\ and k± are the 
momenta parallel and perpendicular to the zone diago- 
nal (magnetic Brillouin zone boundary). The integral of 
Eq. Dl becomes 

1 /-^"""^ de f dk' 

where 9„iin and 9max arc the opening angles of the wedge. 
Thus, the integral over k' diverges at a hot spot. In this 
case, the expansion Eq. C18 must be modified. 
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2. A 

The apparent singularity of A is an artifact. In reality, 
a finite q' shifts the location of the hot spots. To confirm 
this, the hot spot integrals can be evaluated as above. 
First, the shift in the positions of the hot spots is found. 



the points where ^t^ff 



shift by (Eq. 6) kx^a 
in g, the shift is 



k-q/2 



-9/2 ■ 

etc. 



= /i. Denoting the 
then to first order 



2 
1 

2tc 



ilv ^ 
ily 



(D4) 



Expanding about the shifted hot spots (with k' ~ radial 
distance from new hot spots gives 



fit;? 



4:tk'{ag{q) + (3e{q)k'), where ag{q) uq and [3e{q) 
as g ^ 0. Hence, the integral to be evaluated is 



dk'\ 



1 



1 



ot0{q) + f30{q)k' ag + /Sgk' 



(D5) 



which is well behaved as a function of k' and can be 
expanded as a series in q". Note that there are terms 
linear in q, which are cancelled in averaging over the hot 
spots. 

Since there are no singularities, direct numerical eval- 
uation of xil') ~ xW — 0) is straightforward, yieldiing 
A, Fig. 13. [In the numerical calculation, some care is 
needed due to the terms linear in q. While they can- 
cel when summed over all eight hot spots, to accurately 
determine the smaller quadratic term, this summation 
should be carried out at each fc-point prior to summing 
the result for all k.] 



3. A at the C-point 

While the above analysis works for small q near the 
tip of the wedge, it is hard to extend it to the edge of 
the g-plateau, or in particular to the C-point, where the 
plateau width shrinks to zero. It is convenient to intro- 
duce a simplified model*^ , for which the g-dependence of 
X can be calculated analytically. It is convenient to recall 
Fig. 12. While the dashed lines in that figure represented 
an bj shift, they can equally well describe the g-shift of 
the energy denominator, Eq. D2. The plateau edge cor- 
responds to the point where the dashed line intersects 
the Q-shifted FS (horizontal arrows). In the simplified 
model, the energy denominator is linearized, so Ae oc fc^, 
independent of fcy. Chosing q to point along the (7r,7r) 
direction, the FS can be approximated by two circles of 
radius kp, centered at (tt, tt) and (— tt, —it) (for the choice 
of q the other two circles at (tt, — tt) and (— tt, tt) can be 
ignored). The Q-shifted FS is then a circle centered at 
r = (0,0). The FS at (tt, tt) and the Q-shifted FS are 
illustrated in Fig. 29c. 




FIG. 29. (a) Calculated susceptibility x{q) for several val- 
ues of overlap 8. (b) Blowup of plateau region, for Xq ~ Xq^o- 
(c) Model of Fermi surfaces, defining 5, k\ {kj_) and 

Adding the contributions of the overlap of the Q- 
shifted FS with both the FS at (tt, tt) and the one at 



(-TT, -tt), Xq « Iks+q + Ik. 



ks-qj 



ith 



/ = 



''■^ dk I dk 



(D6) 



where the region of integration is over the part of the 
upper FS in Fig. 29c not overlapped by the lower (Q- 
shifted) FS, and k± ranges from zero at the apex of the 
wedge to the middle of the upper FS, kc = kp — ks, where 
ks is the overlap parameter defined in Fig. 29c. To lowest 
order, for ks << kp, 



I = 2^J2kpksyl + l3 + \n\ 



-e(l-/3)(/W3 + ln| 



x/T^-f l' 



(D7) 



with (3 — kc/ks and Q the unit step function. The expres- 
sion for Iks-q must be modified when q > ks and the two 
FSs no longer overlap®"^: Iks-q — 2\/2kp[l — '-f ta,n~^ l/7]i 
with 7 = \/{q — ks)/2kp. The calculated susceptibilities. 
Fig. 29a, display the flat topped plateaus with weak pos- 
itive curvature {A < 0, Fig. 29b). At the plateau edge 
the susceptibility falls sharply, x ~ 1 — 7r7/2 ^ y/q. The 
C-point corresponds to fc^ = 0. 



4. B 

The expression for B may be written exactly as the 
Lu ^ limit of 



B = UReY^i 



1 



/(gfc)-/(gfc+Q) j 

i^k+Q-^k) ii^k+Q-^k)^-'"^' 



(D8) 
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It can be shown that B has a logarithmic correction due 
to the hot spots. The integral can be approximately eval- 
uated by (a) using symmetry to reduce the integral to one 
over an octant of the Brillouin zone containing one hot 
spot, (b) splitting the domain of integration into (i) a 
circle of radius kc about the hot spot, and (ii) the re- 
mainder of the domain, and (c) numerically evaluating 
the integral over domain (ii) while providing an analytic 
approximation to that over (i). Then the k integral over 
the hot spot circle can be written approximately as 



/ = 



(1 - 3P0k)dk 



aik"^ 



1 , 1 



- 

oiekc. 



Oil 



-SPe log 



(D9) 



At T = 0, the integral I must then be integrated in 6 over 

the wedge where the difference in Fermi functions does 
not vanish. The integral from outside the hot spot circle 
will eliminate the fcc-dependence, but should not affect 
the lofj{uj) term. The same logarithmic divergence can 
be found as a byproduct of the calculation of C, below. 

It is difficult to directly evaluate the two-dimensional 
principal value integral for B. Instead, it is much sim- 
pler to evaluate Re{x) via Kramers-Kronig transforma- 
tion of Im{x) a-iid find B by numerical differentiation. 
When this is done, it is found that (a) B is numerically 
very small due to the plateau in iie(x), Fig. lOd, and 
(b) the logarithmic correction is too small to determine 
accurately. 



5. C 

Parameter C is conveniently found by analytically con- 
tinuing Eq. D8 back to the Matsubara frequencies. 



C 



./(efe)-./(efc+Q) 



1 



^^k+Q ^k) ((^ft+d ^k) 



(DIO) 



The wedge integral can be evaluated as for B. The rele- 
vant integral is 



Jo 



dk{l - SPek) 



—tan {—) — m(l 



TT 

2lu„ 



1 

kc 



k 



(Dll) 



{u>„ = uj„/ag) thus giving both the linear in co dissipa- 
tion, and confirming the ln{w) divergence found above. 
Figure 11 shows the divergence of C as /i approaches the 
VHS. Note that it is cut off at increasingly lower frequen- 
cies: the arrows correspond to oJn = ttC. 



Alternatively, Eq. C22 may be used; this can be in- 
tegrated to yield Eq. 14 (in agreement with Sachdev, et 
al.^''), explicitly displaying the divergence at the VHS 
{SxO 0). 



6. u 

Since there is some controvcrsy''^^'^^ concerning u, it 
shall be evaluated in detail. Millis^'* showed that for free 
electrons (parabolic bands) this expression is in general 
well defined, but diverges when Q is a 'spanning' vector 
of the Fermi surface - in the present case, this would 
correspond to the H- and C-points. Abanov, et al.^^ 
found a more severe divergence: u diverges for all /i in 
the hot spot regime. The problem lies in the limit of 
external frequencies 0, momenta or Q. Taking 
this limit on the momenta, the expression for u can be 
written as 



1 



1 



■ (D12) 



The sum over Matsubara frequencies yields 



^ iuiA ^ (ioja — Ae) iiuja 



+ 



iuji {iu)3 — Ae) (iws -|- iiJ2 — Ae) 

K 

1 

{ioji + Ae)(ia;i -|- iui2 + Ae) 

fin+q) . 1 

iL02 {iL03 — Ae){iuji + iui2 + Ae) 
1 

{iuj\ + Ae)(ia;3 -|- iL02 — Ae) 



)], (D13) 



where Ae = eg - e^+Q- Letting uji^± = (wj ± a;i+2)/2 
{i = 1,2), and noting that u)i+ = —0J2+, this simplifies 
to 



(/(e£+Q)-/(ej))VF_ 



where 



W- = -I- Ae). 



(D14) 



(D15) 



Thus in Matsubara frequency space, u is largest for 
wi+ = u)2- = 0, so it should indeed be reasonable to 
estimate it in that limit: 



u{iLOi,0,0) = W 
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,^/(fMl^. (DI6) 

4^ iwi-l-Ae ^ ' 



In turn, it should be possible to approximate u, Eq. D16, 
by its — > limit, if this is nonsingular. Prom Eq. C18, 
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continuation iuji — > 



iCijj + 1 — (5o. Thus, the analytic Es{kz) = —^tgz cosA:2c[cos {k^ + ky)a/2 + cos {kx — ky)a/2] 
j + iS yields = -8t,.^ cos fc^ccos fc:ra/2 cos /cj;a/2, 



w(0,0,0) = C/2 lim 



2BU. 



(D17) 



Due to the plateau in x((5,w), B (Table I) and hence 
u are extremely small. The smallness of u is true only 
in the limit that all external frequencies are small, which 
means that a more complicated expression should be used 
to evaluate u. Moreover, there is an additional problem: 
as found above, B has a correction in In (cj), which would 
formally be divergent. Hence, the model is not fully self- 
consistent, and u will be treated as an empirical param- 
eter. The weak logarithmic divergence will be neglected, 
and u approximated by a constant. 



APPENDIX E: INTERLAYER COUPLING 

1. Dispersion of tz'. Direct and Staggered Staclting 

Andersen, et al.** demonstrated that the anomalous 
form of interlayer hopping in the cuprates, tz = tzo{cx — 
CyY, could be understood by coupling the CU(j2_^2 and 
Op orbitals to the Cn^s orbitals, which have significant 
interlayer coupling. Here, I provide a simplified calcu- 
lation including only these orbitals, and show how the 
dispersion is modified by staggered stacking of the Cu02 
layers. For uniform stacking (Cu above Cu), the hopping 
matrix becomes 
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with Si ~ smkia/2. Here the first (last) row is for the 
Curf2_(i2 (Cu4s) orbital, and the middle rows are for the 
Opx and Opy orbitals, with E^z — — 4t,,2 cos /c^c. In the 



limit A, + Esz 
has dispersion 



>> A >> t,tps, the antibonding band 
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so if tsz << A^, the interlayer hopping has the form 
tzocoskzc{cx - CyY, with tzo = -IQt^tl^t.z/ A? A^. 
While this form had been suggested earlier®^ and found 
experimentally for the bilayer splitting in BSCCO^", it 
should be noted that it is only approximate, and that, 
at least in YBCO, there is considerable splitting of the 
bilayer bands along the zone diagonal**. Nevertheless, 
this form is adequate for the present purposes. 

When successive layers are staggered, the only modi- 
fication to the hopping matrix is in the form of Es{kz), 
which now acquires an in-plane dispersion. 



which leads to Eq. 24. 

2. Estimation of tz from Resistivity Anisotropy 

The dc conductivity can be estimated 
2e2 



(E4) 



i ~ x,y, z, with Q the unit cell volume, Vi — h^^de^/dki, 
and Tj: the scattering rate. Recent ARPES data suggest 
that, when bilayer splitting is resolved, is relatively 
isotropic over the Fermi surface^^. Taking indepen- 
dent of fc, the conductivities are given by integrals over 
the Fermi surface. Figure 30a shows a normalized con- 
ductivity ratio. 



at^ Ozz 



(E5) 



while Fig. 30b shows the resulting normalized interlayer 
hopping izo = tzo\/c/a, which would be required to pro- 
duce a resistivity anisotropy Pzz/Pxx — 1000. For sim- 
plicity, it is assumed that tzo is small, and Ozzl^xx is 
evaluated in the limit i^o ~* 0. It can be seen that (a) 
the staggered stacking reduces the conductivity by ap- 
proximately a factor of 20, independent of doping (except 
near the VHS), so (b) assuming the resistivity anisotropy 
is 1000 for optimally doped LSCO, it is estimated that 
tzidji = 0.11 (for staggered stacking) or 0.025 (for uniform 
stacking). 
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FIG. 30. (a) Normalized conductivity ratio, (Tzz/cr^^a; vs 
doping Sf, for uniform (solid line) and staggered stacking 
(long dashed line and short dashed line, (x20)); and (b) re- 
sulting normalized interlayer hopping t^a for staggered (solid 
line) and uniform stacking (long dashed line and short dashed 
line, (x4.5). 
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3. z-Component of Ordering Vector 

Given a finite interlayer hopping , the first issue is to 
identify the three-dimensional ordering vector: what 
minimizes the free energy? At mean field level, the initial 
magnetic instability will be associated with the state for 
which the RPA denominator first diverges, i.e., the state 
with the largest value of RexoiQ, Qz)- (Note that these 
calculations implicitly assume that the two-dimensional 
ground state involves commensurate order at Q.) For 
uniform stacking, a complicated dependence on doping, 
temperature, and tz is found. Figures 31,32 plot xo vs 
chemical potential for T ~ lOOK, lOK, respectively. The 
shift of the susceptibility peak with doping can readily 
be understood by comparison with Fig. 1. Both temper- 
ature and interlayer coupling act to smear out the VHS, 
and in both cases cause the susceptibility peak to shift to 
smaller chemical potential (lower hole doping). Fig. 3 Id. 
Note that the peak shifts at different rates for different 
Qj-values, showing that the band is developing a consid- 
erable c-axis dispersion. The fastest shift (short dashed 
line in Fig. 31d, corresponding to Qz = 0) can thus be 
considered as representing a crossover from quasi-two- 
dimensional to fully three dimensional dispersion. 



1.4 



1.2 
1.0 





1 1 


1 

(a) ; 


'/ / <-0.10 
/ / i-0.20 
- / -0.30 




'zO , t . " 


'/, 1 


1 0.2 , 0.4 



-0.30 



-0.20 -0.10 




-0.30 -0.20 -0.10 



^(eV) 



FIG. 31. Xo{Q,Qz) at T = WOK vs chemical potential fi, 
for uniform stacking and Qz = n (a), 7r/2 (b), and (c). 
The various curves correspond to tzo/t = 0.01, 0.02, 0.05, 0.1, 
0.2, and 0.5, with the peak in xo shifting to the right with 
increasing tzo- Inset (d): position of peak, fj,max, vs t^o for Qz 
= TT (solid line), n/2 (long dashed line), and (short dashed 
line) . 
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FIG. 32. xo{Q,Qz) vs chemical potential jj., as in Fig. 31, 
but at T = lOK. 

This dispersive shift of the peak in xo leads to a 
doping dependence of the optimal Qz, as illustrated in 
Fig. 33 for tzo = 0.2t. For large hole doping, near the 
tzo ~ VHS, the susceptibility maximum corresponds to 
Qz — tt/c, while near the susceptibility peak, the spin 
modulation becomes incommensurate (intermediate val- 
ues of Qz have the largest susceptibility). There is a 
rapid evolution of the optimal Qz, and beyond the peak 
regime, over essentially the entire electron-doped regime, 
the optimal Qz is 0. This same pattern is repeated for 
smaller i^o , with only the region of the susceptibility peak 
changing. The results are essentially independent of the 
sign of tz ■ 
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FIG. 33. Xo(Q,Qz) vs chemical potential fj., for uniform 
stacking and t^o = 0.2t, and T = lOK (a), or lOOK (b), 
with Qz/tt = 1 (solid line), 0.75 (long dashed line), 0.5 (short 
dashed line), 0.25 (dotted line), (dot-dashed line). 



4. Calculation of 

a. Uniform Stacking 

Given tz and Qz, the parameter Az of Eq. 2 can be 
evaluated: Ux{Q + qzZ,uj = 0) = UxiQ + QzZ,0) + 
A-z{qz — Qz)^ ■ The dominant ordering vectors, Qz — / c 
and Qz = 0, can be analyzed in more detail. For the 
former choice, 
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with fL = - fj:)/kBT, f'J = -fL{l - 2fj:)/kBT, 

Cz — COS kzC, Sz — sinfczC. For the latter case 
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Figure 34 ( 35a) shows how XoiQ, Qz) varies with Qz 
for tzo — O.lt (0.02i), for a number of different dop- 
ings. For the entire electron-doped regime, the peak is at 
Qzm = (Fig. 34b, 35d), crossing over to Qz„i — n/c \n 
the hole doped regime. Away from the peak, the sus- 
ceptibiUty varies as Azq'^, with qz = Qz — Qzm, and 
in the electron-doped regime the full variation can be 
approximated by a cosine. The amplitude of the co- 
sine falls to zero as the C-point is approached. In the 
quasi-two-dimensional regime this amplitude scales with 
i^o- Figure 35b, c shows plots of the best parabolic fit to 
A'^ — Az/c^ for tzo/t = 0.02 (squares) and 0.1 (triangles). 
For tzo/t = 0.1, an alternative A'^ is shown, found by fit- 
ting the full susceptibility as a cosine in qz (circles) . The 
good agreement between the two techniques shows that 
this is a reasonable approximation in the electron-doped 
regime {—0.2eV < /x < 0). Near the susceptibility peak, 
the variation is nonsinusoidal, and the parabolic fit leads 
to a large value for A'^. 
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FIG. 34. {&)xo{Q,Qz) vs Q, for t,o = O.lt, and 
T = lOK, and a variety of chemical potentials fi 
= -0.003559 (solid line), -0.08898 (long dashed line), 
-0.1779 (short dashed line), -0.2669 (dotted line), -0.2847 
(dot-dashed line), -0.2954 (long- long-short-short-short dashed 
line), -0.3025 (long-short-short dashed line), -0.3203 
(dash-dot-dot line), -0.3381 (long-short dashed line), and 
-0.3559 meV (long-short-short-short dashed line). (b) 
= Xo{Q,Qz) — Xo{Q,Qz = 0), where the curves have 
the same meaning as in frame (a). 
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FIG. 35. (a)xo(Q,Q.) vs for t,o = 0.02t, and 
T = lOK, and a variety of chemical potentials fi = -0.003559 
(solid line), -0.08898 (long dashed line), -0.1779 (short 
dashed line), -0.2669 (dotted line), -0.3025 (dot-dashed line), 
-0.3381 (long-long-short-short-short dashed line), -0.3417 
(long-dashed-dotted line), -0.3452 (long-short-short dashed 
line), -0.3488 (long-short-short-short dashed line), and -0.3559 
meV (long-dash-dot-dotted line). (b,c) A'^ = A;^/Uc^ vs /i for 
tzo/t = 0.02 (squares, A'^ x 25) and 0.1 (triangles,circles) . (d) 
Qzm VS fj, for t^o/t — 0.02 (squares) and 0.1 (triangles). 



FIG. 36. {si)A'^ — Az/U(? vs chemical potential ^ for Qz ~ 
(solid lines) or tt (dashed lines) , for a variety of values of t^o 
and T = IQQK. In order of increasing amplitude, the values 
are tzo/t = 0.01, 0.02, 0.05, 0.1, 0.2, and 0.5. (b) Scaling 
of A'j°^ with {tzo/tf. Curves are tzo/t = 0.01 (solid line), 
0.02 (long dashed line), 0.05 (short dashed line), 0.1 (dotted 
line), 0.2 (dot-dashed line), and 0.5 (dot-dot-dashed line), (c) 
Comparison of max{Az) for staggered stacking (solid line) 
and uniform stacking (triangles, xl/20) at tzo/t ~ 0.1. 



b. Staggered Stacktng 

The same calculations can be repeated for the of 
Eq. 24, associated with staggered stacking; Fig. 36a 
shows Az calculated from Eqs. E6, E7 at = (solid 
lines) and tt (dashed lines). The frustration induced by 
staggering of the Cu02 layers is reflected in a strong sup- 
pression of the g^-dependence of x, which leaves a small 
residual contribution quadratic in tzo, Fig. 36b. Since tz 
vanishes at (tt, 0), there is no shift of the susceptibility 
peak with doping. Note the symmetry of the Az values 
between and tt. In fact, xiQz) is closely sinusoidal, 
particularly for small tzo, with maxima either at tt or 
0. Thus, near either the H- or C-points, the maximum 
of X corresponds to Qz — tt. For intermediate dopings, 
Qz = is favored. At two distinct chemical potentials, 
the amplitude of the cosine collapses and changes sign. 
At the crossing points, x is independent of Qz, leading 
formally to Tn 0. Note from Fig. 36c that the sup- 
pression of Az is approximately in the same ratio as that 
of the resistivity, found above. 
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5. Calculation of Tjv 

When there is a finite interlayer hopping tz 
becomes 



Eq.C37 
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where 



5 + Az_q'z^ and T3D = 7r^Az/2Ce'^c^. (A 
small correction to 60 is neglected. Treating the Qz de- 
pendence as a cosine rather than a cutoff quadratic leads 
to qualitatively similar results.) Thus a finite Az always 
cuts off the divergence found in Eq.C37, leading to a fi- 
nite Tn whenever there is a zero-temperature Neel state 
(e.g., up to a QCP). It should be noted that the above 
calculation implicitly assumed that T > T^r) ~ Az'. for 
T < T^D the logarithm is cut off and the system be- 
haves like an anisotropic three-dimensional magnet. For 
tzo/t < 0.1, the system is generally in the quasi- two- 
dimensional limit, Fig. 37a. Figure 37b compares the 
mean-field Neel transition with the Neel transition found 
assuming uniform stacking and finite interlayer couplings 
tzo/t = 0.1, 0.02, and 2 x 10"^ [the last found by scahng 
the T^D for tzo/t — 0.02 by the ratio of tl^s]. It is seen 
that T/v ^ as ^ 0, albeit exceedingly slowly. 
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FIG. 37. {a)T-iD vs fi for T = lOK and uniform stack- 
ing with tzo ~ O.lt (circles, xl/25) or 0.02 (squares), or 
staggered stacking with f20=0.02 (triangles), (b) Tn vs x, 
comparing mean field transition (solid line) with interlayer 
coupling models (uniform stacking) assuming tzo/t = 0.1 
(long dashed line), 0.02 (short dashed line), and 2x10"® 
(dot-dashed line) , and the staggered stacking model assuming 
tzo/t ~ 0.1 (dot-dot-dash line). 



FIG. 38. (a) Mean field Tjv vs x assuming paramagnetic 
Ueff (Appendix B). (b) Corresponding Tm vs x, calculated 
using Eq. E8. Squares — staggered stacking with tzo/t = 0.1; 
triangles — uniform stacking with tzo/t — 0.02; solid line and 
circles = data, as in Fig. 25. 



The above calculations are for uniform stacking. For 
staggered stacking is reduced, in approximately the 
same ratio as the resistivities. Hence, the staggered 
stacking with tzo/t = 0.1 should be comparable to uni- 
form stacking with tzo/t = 0.02, as observed. Fig. 37. 
While Tn technically goes to zero for staggered stacking 
near x = —0.0838, the decrease is logarithmic, and in 
practice no more than a weak dip is expected to be ob- 
served (the point with T/v = OK is omitted from the plot 
in Fig. 25). Hence, if tzo is estimated from the resistivity, 
it will be nearly impossible to distinguish uniform from 
staggered stacking via measurements of Tjv. 

In the above calculations, a constant value of A was 
assumed for each doping, as given in Fig. 13. In fact, 
for the electron-doped cuprates, A ^ 1/T^-^ for T > T4, 
Fig. 9. This would cause an enhancement of the logarith- 
mic correction, ~ T^-^, tending to pin Tn close to TjJ. 
For the present parameter values, this could reduce Tpf 
by roughly a factor of two, still larger than the experi- 
mental values. 

A more likely source of the discrepancy is the possible 
temperature dependence of Ueff, Appendix B. The large 
Ueff at half filling arises from lack of screening, in the 
presence of a Mott gap - and is appropriate in analyzing 
the low-T Fermi surfaces found in ARPES. For calculat- 
ing the onset of the Mott gap, the mean field T^, it is 
more appropriate to use the paramagnetic susceptibility, 
as in Fig 28a. When this is done, considerably smaller 
transition temperatures are found, both at the mean field 
level. Fig. 38a, and when fluctuations and interlayer hop- 
ping are included. Fig. 38b. While the latter are closer 
to the experimental values, no attempt has been made to 
correct Ueff for the short range gap. 
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